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EVALUATION 


1.  As  future  tactical  and  strategic  sensors  become  sore  automated  and  unattended, 
their  ptocessing  will  increasingly  treat  the  environmental  scatter  as  information 
rather  than  "clutter"  or  interference.  By  inference  from  measureable  quantities 
and  statistics,  the  processor  will  recognize  the  existence  of  weather,  diaff,  dis¬ 
crete  targets ,  homogeneous  distributed  areas ,  shadowing  as  opposed  to  specular  re¬ 
flection,  and  other  environmental  categories.  This  information  will  allow  the  sys¬ 
tem  to  adapt  its  waveform,  energy  budget,  and  detection/CFAR  algorithms  for  optimum 
performance.  This  effort  investigated  algorithms  which  may  be  usable  in  such  a 
recognition  scheme. 

2.  Stress  was  placed  on  the  strengths,  weaknesses,  and  peculiarities  of  two  prom¬ 
ising  algorithms  when  applied  to  the  recognition  of  shadowed  areas.  Ihe  program 
has  provided  an  excellent  insight  into  the  robustness  and  sensitivities  to  data 
format  of  each  algorithm  and  the  care  required  in  the  interpretation  and  use  of 
algorithmic  results.  This  report  will  be  a  useful  reference  in  the  development 

of  clutter  model  and  associated  adaptive  processing  for  future  intelligent  unat¬ 
tended  systems. 

UJjyt** j 

WILLIAM  L.  SIMKINS 
Project  Engineer 


MODELING  OF  SHADOWS  IN  RADAR  CLUTTER 

1.0  INTRODUCTION 

This  report  presents  the  results  of  the  work  done  during 
the  period  June-December  1979  on  investigating  the 
applicability  of  nonparametric  clustering  techniques  to  the 
extraction  of  features  from  radar  clutter.  The  Forgy-Jancey- 
[1]  and  minimum  spanning  tree  [2]  algorithms  were  investigated, 
with  particular  application  to  the  problem  of  locating  and 
characterizing  shadowed  regions  in  ground  clutter.  The 
significant  result  of  the  work  is  an  exposition  of  the 
strengths  and  weaknesses  of  the  two  clustering  techniques  when 
applied  to  both  simulated  and  actual  radar  returns.  It  is 
expected  that  the  same  strengths  and  weaknesses  should  appear 
in  applications  involving  extraction  of  clutter  features  other 
than  shadows. 

The  motivation  for  undertaking  the  present  study  of  ground 
clutter  is  summarized  in  [33.  The  need  has  existed  since 
before  the  Second  World  War  to  extend  the  ability  of  radars  to 
detect  and  identify  targets  of  low  cross-section  in  severe 
clutter.  Initial  efforts  were  concerned  with  increasing 
transmitter  power,  but  for  the  past  two  decades  the  emphasis 
has  been  on  development  of  the  signal  processing  abilities  of 
radar  receivers  [4].  The  recent  explosive  developments  in 
digital  technology  are  expected  to  result  in  a  new  generation 
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of  sophisticated  signal-processing  radars,  able  to  take 
advantage  of  random-access  memories  and  high-speed 
computational  ability.  The  need  to  specify  how  such  radars  are 
to  be  designed  and  used  provides  a  requirement  for  an  improved 
understanding  of  the.  structure  of  backscatter  clutter. 

There  are  several  indications  that  the  present 

understanding  of  backscatter  clutter  is  inadequate.  First, 
clutter  measurement  efforts  concerned  with  various  terrain  and 
cultural  features  have  produced  data  of  only  limited 

consistency.  In  particular,  it  is  not  uncommon  to  find  a  15  to 
20  dB  spread  in  measurements  of  intensity  returns  made  at 
different  times  from  supposedly  identical  terrain.  Second, 
there  has  been  only  limited  success  in  predicting  backscatter 
from  toptgraphical  maps,  even  though  the  terrain  relief, 
surface  features,  and  ground  cover  are  specified  on  the  maps. 
Finally,  radars  employing  signal  processing  systems,  such  as 
CFAR,  that  are  based  on  particular  clutter  models  [53  have 
frequently  been  found  to  perform  more  poorly  than  analysis 
based  on  the  clutter  models  would  suggest. 

There  are  several  possible  reasons  for  the  poor 
consistency  of  measurement  data  and  the  difficulty  of 
predicting  backscatter  intensity  from  terrain  descriptions. 
First,  it  is  likely  that  backscatter  is  highly  sensitive  to 
surface  parameters  such  as  moisture,  temperature,  and  wind 
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velocity,  that  are  difficult  to  measure  simultaneously  with  the 
radar  return,  and  have  not  always  been  included  as  a  part  of 
the  recorded  "ground  truth".  Second  is  the  possibility  that 
the  clutter  returns  vary  with  the  parameters  of  the  radar 
itself  in  ways  that  have  not  yet  been  clearly  recognized. 

Much  of  the  radar  signal  processing  effort  to  date  has 
been  based  on  the  assumption  that  clutter  intensity  can  be 
modeled  by  a  random  variable.  The  motivation  to  develop  such  a 
model  is  strong,  because  of  the  success  of  design  methods  using 
the  Neyman-Pearson  decision  rule  in  connection  with  a  gaussian 
or  rayleigh  model  for  thermal  noise.  Such  decision-theoretic 
methods  could  be  easily  extended  to  detection  of  targets  in 
clutter,  if  only  the  probability  density  function  of  the 
clutter  intensity  could  be  found.  It  is  likely,  however,  that 
ground  clutter  is  too  heterogeneous  to  be  represented  by  a 
single  random  variable,  or  possibly  even  by  a  single  type  of 
random  variable.  It  is  therefore  necessary  that  an  examination 
of  the  structure  of  radar  clutter  be  carried  out  in  a  way  that 
does  not  rely  on  knowledge  of  a  probability  distribution. 

Methods  for  investigating  the  structure  of  unknown  data 
have  been  developed  in  connection  with  research  in  taxonomy 
[1,6]  and  chemistry  [73.  These  methods  are  referred  to  as 
either  classification  or  slU5.tfirlng  methods,  depending  on 
whether  or  not  a  statistical  structure  is  assumed  for  the 
unknown  data.  The  classification  methods  assume  that  each  data 
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entity  belongs  to  one  of  several  possible  categories,  each 
category  characterized  by  a  probability  density  function.  Most 
classification  algorithms  require  training  on  data  whose 
classification  by  category  is  known.  The  training  allows  the 
algorithm  to  estimate  the  density  functions  corresponding  to 
the  various  categories.  The  estimated  density  functions  are 
then  used  in  classifying  future  unknown  data.  In  contrast, 
clustering  methods  do  not  assume  that  categories  of  data  can  he 
characterized  by  density  functions.  Instead,  data  entities  are 
grouped  into  clusters  by  examining  the  distances  between 
entities  based  on  some  suitable  measure  of  distance.  A 
detailed  discussion  and  comparison  of  many  available 
classification  and  clustering  methods  is  presented  in  [6]. 

In  the  present  research  the  view  has  been  taken  that 
ground  clutter  from  arbitrary  sources  can  not  be  reliably 
represented  in  terms  of  probability  density  functions  [8], 
Attention  was  therefore  focussed  on  an  investigation  of  the 
various  clustering  methods.  The  research  has  proceeded  along 
three  paths:  first  an  investigation  of  the  Forgy-Jancey 

clustering  algorithm,  second  an  investigation  of  the  minimum 
spanning  tree  algorithm  contained  in  the  program  package 
ARTHUR,  and,  finally,  an  investigation  into  some  pre-processing 
techniques  applicable  to  either  of  the  clustering  algorithms. 
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A  description  of  the  Forgy-Jancey  algorithm  along  with  the 
results  obtained  from  applying  it  to  a  variety  of  simulated  arid 
actual  clutter  returns  is  presented  in  Section  2.  Although  the 
Forgy-Jancey  algorithm  is  well-known,  it  was  necessary  to 
create  our  own  programs,  and  to  devise  from  scratch  the 
necessary  input-output  and  display  routines.  Investigation  of 
the  minimum  spanning  tree  clustering  method  is  described  in¬ 
detail  in  Section  3.  This  investigation  began  with  acquisition 
of  the  ARTHUR  program  package  [7]  from  the  University  of 
Georgia.  The  necessary  job-control  language  was  created,  and 
the  ARTHUR  package  compiled  so  that  it  is  currently  available 
on  the  Clarkson  College  computer.  The  minimum  spanning  tree 
algorithm  was  tested  on  some  of  the  same  simulated  data  as  was 
devised  for  the  Forgy-Jancey  tests.  Section  4  contains  the 
results  of  some  pre-processing  efforts  that  were  tried  on 
actual  clutter  data.  The  pre-processing  was  motivated  by  the 
need  to  reduce  the  quantity  of  data  to  a  managable  amount  for 
the  ARTHUR  package.  The  pre-processing  is  typical  of  the  steps 
that  would  be  taken  in  an  actual  clustering  effort,  and 
provided  an  indication  of  the  amounts  of  CPU  time,  programming 
time,  and  costs  that  are  involved  in  such  steps.  Section  5 
contains  an  overall  summary  and  the  conclusions. 

Sections  1,  2,  and  5  were  written  by  Bruce  A.  Black, 
Section  3  by  William  Ladew,  and  Section  4  by  Mohammed 
Arozullah. 


5 


MODELING  OF  SHADOWS  IN  RADAR  CLUTTER 


2.0  THE  FORGY-JANCEY  ALGORITHM 


2. 1  The  Algorithm 

This  section  discusses  the  operation  of  the  Forgy-Jancey 
clustering  algorithm,  and  the  tests  that  were  made  on  simulated 
and  actual  radar  data  to  probe  the  strengths  and  weaknesses  of 
the  algorithm  as  a  tool  for  identifying  the  structure  of  ground 
clutter  from  shadowed  terrain.  The  computer  program  for  the 
Forgy-Jancey  algorithm  was  obtained  from  Fordon  [91.  Several 
modifications  to  the  algorithm  were  made  to  allow  greater 
flexibility  in  manipulating,  saving,  and  displaying  the  input 
and  output  data.  The  algorithm  employed  in  all  of  the  tests  is 
due  to  Forgy  [10],  Jancey's  version  [1]  contains  a  heuristic 
to  avoid  local  minima,  but  was  not  needed  in  the  present 
application.  The  basic  Forgy-Jancey  program  along  with  the 
modifications  is  presented  in  block-diagram  form  below.  A 
listing  of  the  important  subroutines  is  presented  in  an 
Appendix. 

Figure  2.1  is  a  flowchart  showing  the  operation  of  the 
Forgy-Jancey  algorithm.  This  algorithm  forms  clusters  of  data 
entities,  where  each  data  entity  is  a  vector  of  NV  components. 
The  number  NV  of  components  in  a  data  entity  (number  of 
variables),  the  number  NE  of  data  entities,  the  number  NC  of 
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clusters,  and  the  coordinates  of  the  cluster  centroids  must  all 
be  specified  by  the  user.  In  this  version  of  the  algorithm  the 
number  NC  of  clusters  cannot  be  modified  by  the  program;  the 
program  can  neither  create  nor  delete  clusters.  In  operation, 


the  program 

will  partition 

the  NE 

data 

entities 

into 

NC 

clusters.  A 

measure  of  the  compactness 

of  the 

clusters 

that 

is 

used  as  a 

figure  of 

merit  by 

the 

program 

is 

the 

total  summed 

deviations.  For 

a  given 

partition  of 

the 

data 

into  clusters,  this  is  the  sum  over  all  the  entities  of  the 
distance  between  each  entity  and  the  centroid  of  the  cluster  to 
which  it  is  assigned. 

Forgy's  algorithm  begins  with  an  initial  set  of  NC  cluster 
centroids.  In  the  main  loop  of  the  program,  each  data  entity 
is  assigned  in  turn  to  the  nearest  cluster  centroid.  When  all 
of  the  entities  have  been  assigned,  the  cluster  centroids  are 
recomputed.  Each  iteration  requires  a  complete  pass  through 
the  data.  It  is  easily  shown  that  a  pass  through  the  main  loop 
cannot  increase  the  value  of  the  total  summed  deviations,  so 
that  the  algorithm  must  converge  to  at  least  a  local  minimum. 
In  each  pass  through  the  main  loop  the  program  counts  the 
number  of  entities  whose  cluster  assignment  is  changed.  When 
the  number  of  changes  falls  below  a  user-supplied  parameter 
MINREL,  the  program  halts.  Except  where  otherwise  noted, 
programs  were  run  with  MINREL  equal  to  zero,  corresponding  to 
complete  convergence. 
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The  distances  between  entities  and  cluster  centroids 
required  by  the  algorithm  are  computed  by  a  user-supplied 
subroutine  DIST.  Because  this  routine  is  externally  supplied, 
the  distance  function  can  be  configured  to  meet  the  specific 
application.  In  particular,  it  is  possible  to  base  the 
distance  calculation  on  a  subset  of  the  components  of  the  data 
entities.  Both  euclidean  and  non-euclidean  distance  measures 
are  allowed. 

The  program  path  for  the  Forgy-Jancey  algorithm  is  shown 
in  Figure  2.2.  The  program  MAIN  and  the  subroutines  EXEC, 
USER,  RESULT,  KMEAN,  and  USROUT  were  all  supplied  or  modified 
as  a  part  of  the  present  effort.  The  program  supplies  as 
outputs: 

1.  A  raw  membership  list,  listing  cluster  membership  for 
each  data  entity. 

2.  A  sorted  membership  list,  giving  for  each  cluster  the 
sequence  numbers  of  its  members. 

3.  A  sorted  data  file.  The  data  entities  themselves  are 
sorted  into  groupings  by  cluster  membership.  This  file 
can  be  saved  on  disk,  thus  allowing  the  entities  forming  a 
single  cluster  to  be  isolated  for  further  processing. 

4.  The  program  prints  the  final  centroid  locations  for 
each  cluster,  the  number  of  data  entities  belonging  to 
each  cluster,  and  the  value  of  the  total  summed 
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2.2  Testing 

For  initial  testing  of  the  Forgy-Jancey  algorithm,  a  set 
of  10,000  data  entities  representing  a  100-point  by  100-point 
"scene”  was  generated.  Each  data  entity  was  a  three-vector,. 
(x,y,I),  where  x  and  y  represent  the  position  of  the  entity  in. 
rectangular  coordinates  and  I  represents  the  intensity  of  the- 
radar  return  from  position  (x,y).  Two  rectangular  regions  were 
set  aside  in  the  100x100  array  to  represent  "shadows",  one 
region  occupying  651  points  and  the  other  occupying  121 
points.  The  test  scene  is  depicted  in  Figure  2.3. 

For  the  initial  test  of  the  algorithm,  the  value  0  was 
given  to  the  intensity  component  of  data  entities  in  the 
shadowed  areas,  while  the  value  1  was  given  to  the  intensity 
component  of  data  entities  in  the  unshadowed  area.  The 
algorithm  was  applied  with  the  number  of  clusters  NC  set  equal 
to  two.  The  distance  function  DIST  was  set  to  compute  the 
distance  between  two  entities  as  the  magnitude  of  the 
difference  between  their  intensities.  The  output  of  the 
program  is  shown  in  Figure  2.4.  The  points  marked  1  correspond 
to  cluster  1,  while  the  points  marked  with  a  dot  correspond  to 
cluster  2.  Note  that  cluster  1  contains  exactly  the  entities 
in  the  shadowed  array, 

A  second  clustering  was  carried  out  using  only  the 
entities  in  cluster  1  as  data.  This  time  the  function  DIST  was 
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changed  to  compute  the  distance  between  two  entities  as  the 
euclidean  distance  between  their  positions.  The  program  was 
run  with  NC  =  2.  The  output  is  shown  in  Figure  2.5.  In  this 
figure  the  points  marked  with  a  dot  were  not  included  in  the 
input  data.  The  points  marked  1  and  2  represent  the  two 
clusters  obtained  by  the  algorithm.  Note  that  these  correspond 
exactly  with  the  two  shadows. 

A  second,  more  demanding,  test  of  the  Forgy-Jancey 
algorithm  was  made  using  random  numbers  as  simulated  intensity 
values.  A  program  was  written  to  produce  a  sequence  of  numbers 
having  a  rayleigh  distribution.  The  distribution  parameter 
sigma  is  entered  as  an  input  tr>  the  program.  To  establish  the 
test  data,  10,000  rayleigh-distributed  numbers  were  generated 
with  sigma  s  1.  These  numbers  were  assigned  as  intensity 
values  to  the  data  entities  in  the  100  by  100  array.  Next  a 
sequence  of  rayleigh-distributed  numbers  was  generated  with 
sigma  much  greater  than  one.  These  numbers  were  added  to  the 
intensity  values  of  all  the  data  entities  outside  the  two 
shadowed  areas.  Thus  a  set  of  test  data  was  established  to 
represent  two  shadows  containing  only  rayleigh-distributed 
thermal  noise  in  a  larger  region  of  brighter 
rayleigh-distributed  clutter. 

As  in  the  previous  test,  the  Forgy-Jancey  algorithm  was 
applied  twice,  once  to  separate  the  low-intensity  shadows  from 
the  bright  clutter  background,  and  a  second  time  to  separate 
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» 

i 

the  low-intensity  cluster  into  individual  shadows.  The  result  < 

the  first  clustering  is  shown  in  Figure  2.6.  The  bright  j 

clutter  in  this  case  has  sigma  set  to  10,000.  In  the  figure,  1 

the  dark  cluster  is  indicated  by  ones  and  the  bright  cluster  by  j 

t 

twos.  Note  that  the  dark  cluster  includes  many  points* 
belonging  to  the  bright  region.  In  fact,  the  algorithm  was 
unable  in  this  case  to  locate  the  dark  region;  the  program  was 
instructed  to  find  two  clusters,  and  it  did  so  by  splitting  the 
bright  data  half.  The  shadow  data  are  simply  included  in  the 

j 

darker  of  the  two  clusters.  I 

i 

To  provide  a  greater  apparent  separation  between  the  dark  \ 

I 

and  bright  data,  and  to  provide  an  example  of  the  use  of  | 

i 

coordinate  scaling,  the  program  was  modified  to  interpret  { 

intensity  values  in  decibels.  This  was  effected  by  employing  j 

the  program  USER  to  replace  the  intensity  values  I  for  each  ] 

entity  by  their  logarithms  lnl  as  the  entities  were  read  into 
core.  The  result  of  a  clustering  done  on  the  rayleigh  data 
using  this  logarithmic  transformation  is  shown  in  Figure  2.7. 

Perfect  separation  of  the  dark  and  bright  regions  was 
obtained.  Figure  2.8  was  obtained  using  only  the  dark  cluster 
as  data.  The  dark  cluster  was  split  on  the  basis  of  position 
into  the  two  clusters  shown  by  points  marked  1  and  2. 

The  gap  between  the  dark  and  bright  intensity  values  can 
be  reduced  by  changing  the  value  of  sigma  used  in  generating 
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the  bright  rayieigh  data.  To  create  a  more  demanding  exercise 
for  the  algorithm,  the  data  set  was  recreated  with  sigma  equal 
to  50.  Logarithmic  scaling  of  the  intensity  values  was  again 


used. 

Figure  2.9 

shows 

the  result  of 

clustering 

for  light 

and 

dark 

intensity  regions. 

Points  in 

the 

dark 

cluster 

are 

indicated  by  ones, 

and 

points  in 

the 

bright 

cluster 

are 

indicated  by  twos. 

Note 

that  a  number 

of 

bright 

points 

have 

been 

misclassified 

by 

being  included 

in 

the  dark 

cluster. 

It 

would 

be  expected 

that 

statistical 

variation 

among 

the 

intensity  values  of  the  non-shadow  data  would  lead  to  some 
"bright"  points  with  intensity  values  that  are  in  fact  as  low 
as  the  values  of  the  intensities  of  the  shadow  data.  The 
number  of  misclassif ied  points  shown  in  the  figure  corresponds 
roughly  with  the  number  that  would  be  expected  on  statistical 
grounds. 

Problems  with  convergence  of  the  Forgy-Jancey  algorithm 
were  encountered  for  the  first  time  in  generating  the  clusters 
shown  in  Figure  2.9.  To  generate  this  figure  the  algorithm  was 
halted  using  MINREL  s  100;  i.e.  on  the  first  pass  for  which 
fewer  than  100  data  entities  were  reassigned.  Figure  2.10 
shows  the  result  of  running  the  algorithm  all  the  way  to 
convergence,  a  process  which  required  nearly  five  minutes  of 
CPU  time.  It  is  interesting,  and  somewhat  surprising,  to  note 
that  the  dark  cluster  of  Figure  2.10  includes  even  more 
misclassified  "bright”  points  than  does  the  dark  cluster  of 
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Figure  2.9. 

The  dark  cluster  of  Figure  2.10  was  used  as  the  data  base 
for  the  second  clustering  pass,  in  which  individual  shadows 
were  separated.  Figure  2.11  shows  the  result  of  clustering  by 
position  with  NC  equal  to  two.  In  this  case  the  two  shadows, 
were  correctly  separated,  but  each  cluster  consists  of  a  dense- 
core,  corresponding  to  the  actual  shadow,  surrounded  by  a 
sparse  "halo"  of  points.  In  an  attempt  to  separate  the  shadows 
from  their  halos,  several  clustering  passes  were  made  on  the 
data  with  NC  greater  than  two.  Figure  2.12  shows  the  result  of 
clustering  by  position  when  NC  equals  three.  The  points  marked 
zero  are  excluded  from  the  data  set.  The  three  clusters  are 
indicated  by  points  marked  1,  2,  and  3,  respectively.  (Cluster 
two  is  in  the  upper  right-hand  corner  of  the  figure.)  Note  that 
the  actual  shadows  still  form  the  cores  of  separate  clusters, 
but  that  both  are  still  surrounded  by  halos.  Figures  2.13, 
2.14,  and  2.15  are  the  results  of  clustering  with  NC 
respectively  equal  to  four,  five,  and  six.  Figure  2.16  shows 
the  case  with  NC  equal  to  seven.  One  of  the  shadows  has  been 
split  into  two  clusters,  marked  1  and  7  in  the  figure.  This  is 
an  undesirable  result,  from  which  two  lessons  can  be  learned: 
first  that  the  Forgy-Jancey  algorithm  will  find  as  many 
clusters  as  it  is  instructed  to,  whether  or  not  these  clusters 
are  meaningful,  and  second,  that  increasing  the  number  of 


MODELING  OF  SHADOWS  IN  RADAR  CLUTTER 

clusters  is  not  an  effective  method  for  removing  the  halo  from 
a  dense  core  of  points. 

Identifying  the  value  of  NC  that  corresponds  to  the  number 
of  clusters  actually  present  is  difficult  when  dealing  with 
unknown  data,  especially  data  of  high  dimensionality  that 
cannot  be  readily  plotted.  Jancey  [1]  suggests  beginning  with 
NC  r  1  and  running  the  algorithm  repeatedly  with  HC  =  2,  3, 
etc.  As  the  number  of  clusters  increases,  the  value  of  the 
total  summed  deviations  obtained  at  convergence  will  decrease. 
As  "natural"  clusterings  are  obtained,  there  should  be 
significant  drops  in  the  value  of  the  total  summed  deviations. 
An  example  of  this  effect  is  presented  in  [1],  Figure  2.17 
shows  the  values  obtained  for  total  summed  deviations  in  the 
above  rayleigh  example  when  NC  =  1  through  7.  Note  that  the 
most  significant  drop  occurs  when  NC  is  raised  from  one  to 
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2.3  Removal  of  Outliers 

Fordon  [93  suggests  a  modification  of  the  Forgy-Jancey 

algorithm  to  allow  rejection  of  data  entities  that  are  too  far 

from  any  of  the  cluster  centroids.  A  data  entity  can  be 

isolated  if  it  satisfies  the  relation 

}D  -*D  {  <  (OUTLYR)xD  , 

avg  min  avg 

where  D  is  the  distance  between  the  entity  and  the  nearest 
min 

cluster  centroid  and  D  is  the  average  of  the  distances 

avg 

between  the  data  entity  and  each  of  the  cluster  centroids.  The 
parameter  OUTLYR  is  assigned  a  value  between  zero  and  one  by 
the  programmer.  A  value  of  zero  results  in  no  entities  being 
isolated,  while  a  value  of  one  results  in  every  entity  being 
isolated.  The  main  loop  of  the  Forgy-Jancey  algorithm  was 
modified  so  that  whenever  an  isolated  entity  is  encountered  it 
is  designated  as  the  centroid  of  a  new  cluster.  This  allows 
the  program  to  increase  the  number  of  clusters.  A  provision 
was  included  to  limit  the  number  of  clusters  to  a 
programmer-specified  maximum,  since  CPU  time  was  found  to 
increase  dramatically  with  the  number  of  clusters.  In 
interpreting  the  results  of  a  clustering  run,  the  programmer 
has  the  option  of  disregarding  all  newly-formed  clusters  as 
representing  outliers;  i.e.  sets  of  anomalous  entities  that  lie 
outside  the  known  clusters,  or  of  accepting  all  of  the  clusters 
in  the  expectation  that  the  program  has  found  a  more  natural 
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number  of  clusters  than  the  programmer  originally  specified. 

An  example  of  a  clustering  run  using  the  modified 

Forgy-Jancey  algorithm  is  shown  in  Figure  2.18.  The  data  set 
is  the  lOOxIOO-point  array  depicted  in  Figure  2.3,  with  each 
entity  having  an  intensity  value  of  zero  or  one.  Three 
outliers  were  created,  having  intensity  values  of  zero,  and 
lying  outside  the  simulated  shadows.  The  Forgy-Jancey 

algorithm  was  run  twice,  once  to  separate  out  the  low-intensity 
points  and  once  to  cluster  these  points  into  individual 
shadows.  Figure  2.18  is  the  result  of  the  second  run.  To 

generate  this  figure,  the  program  was  started  with  NC  =  2  and 

OUTLYR  s  0.6.  The  program  responded  by  creating  two  additional 
clusters,  identified  in  the  figure  by  the  numbers  3  and  4. 
Although  the  program  successfully  located  the  three  outlying 
points,  it  should  be  noted  that  the  value  used  for  the 
parameter  OUTLYR  is  fairly  critical.  The  program  was  unable  to 
locate  all  of  the  outliers  when  run  with  OUTLYR  s  0.4. 


17 


MODELING  OF  SHADOWS  IN  RADAR  CLUTTER 


2.4  Testing  of  Actual  Radar  Data 

A  tape  containing  a  segment  of  a  PPI  display  was  provided 
by  the  Signal  Processing  Laboratory  at  RADC,  courtesy  of 
William  L.  Simkins,  Jr.  The  tape  contains  data  in  essentially 
the  same  format  as  the  simulated  data  described  above:  each- 


data  entity  is 

comprised 

of 

a  set  of 

three 

components 

representing  x  and 

y  position 

and 

intensity . 

The 

data  were 

created  by  quantizing  a  photograph  of  the  original  CRT  display, 
and  in  consequence  the  intensity  scale  of  the  data  is 

compressed  by  the  responses  of  the  CRT  phosphor  and  the  film. 

* 

Figure  2.19  depicts  one  quadrant  of  the  radar  data.  The  grey 
scale  used  in  the  figure  consists  of  the  symbols 

(W, », A, I, blank) ,  in  order  from  highest  to  lowest 
intensity.  Figure  2.20  shows  the  low-intensity  regions  of  the 
quadrant  of  data.  The  set  of  low-intensity  entities  was  used 
as  a  data  base  for  the  Forgy-Jancey  clustering  algorithm. 
Numerous  runs  were  made,  using  various  values  of  input 

parameters.  A  typical  output  is  shown  in  Figure  2.21.  To 
generate  this  output  twenty-one  initial  cluster  centroids  were 
chosen,  corresponding  roughly  to  the  centers  of  the  clusters 
visible  in  Figure  2.20.  The  outlier  parameter  was  set  to 

OUTLYR  =  0.9,  and  the  algorithm  was  stopped  when  fewer  than  33 

entities  changed  cluster  membership  (MINREL  s  33).  The  program 
ran  for  nearly  four  minutes  and  generated  a  total  of 


MODELING  OF  SHADOWS  IN  RADAR  CLUTTER 


twenty-five  clusters. 

Examination  of  Figure  2.21  reveals  the  twenty-five 
clusters,  designated  by  the  numbers  zero  through  nine  and  the 
letters  A  through  0.  A  close  examination  of  this  Figure 
reveals  that  many  of  the  clusters  found  by  the  program  do  not 
correspond  to  the  natural  clusters  evident  in  Figure  2.20.  In 
particular,  attention  is  called  to  cluster  K,  which  is  composed 
of  the  adjacent  parts  of  four  natural  clusters,  clusters  C  and 
0,  which  split  a  single  natural  cluster,  cluster  2,  which 
contains  some  entities  that  should  properly  belong  to  cluster 
3,  and  cluster  D,  which  encroaches  on  the  territory  of  clusters 
C  and  0. 

In  processing  the  actual  radar  data,  the  Forgy-Jancey 
algorithm  had  a  great  deal  of  difficulty  in  locating  the 
natural  clusters.  It  is  likely  that  this  difficulty  stems  from 
two  sources:  the  highly  irregular  shapes  of  the  natural 
clusters,  and  the  fact  that  the  natural  clusters  vary  greatly 
according  to  size.  It  is  believed  that  the 
total-summed-deviations  figure  of  merit  used  by  the  algorithm 
tends  to  emphasize  the  compactness  of  clusters.  Thus  whenever 
a  cluster  becomes  too  long  and  narrow,  the  algorithm  will  split 
it  into  two  smaller  clusters.  The  procedure  used  by  the 
program  for  generating  new  clusters  tends  to  result  in  clusters 
of  about  equal  size.  Hence  the  difficulty  in  dealing  with  a 
data  set  that  contains  clusters  as  small  as  cluster  7  and  as 
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larGe  as  cluster  6. 
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2.5  Conclusions 

The  Forgy-Jancey  clustering  algorithm  is  a  powerful  tool 
whose  principle  advantages  seem  to  be  its  generality  and 
simplicity.  The  generality  results  from  the  ability  of  the 
algorithm  to  handle  data  sets  of  essentially  any  size, 
consisting  of  data  entities  of  any  dimension.  The  number  of 
clusters  that  the  algorithm  can  handle  is  also  essentially 
unlimited,  and  the  programmer  is  free  to  prescribe  arbitrary 
coordinate  scaling  and  an  arbitrary  distance  measure.  The 
simplicity  of  the  algorithm  is  helpful  when  modifications  are 
needed.  In  carrying  out  the  present  study  it  was  necessary  to 
add  several  input  and  output  routines,  and  to  modify  the  main 
loop  of  the  algorithm  to  incorporate  the  ability  to  generate 
new  clusters. 

The  drawbacks  to  using  the  Forgy-Jancey  algorithm  are 
clearly  illustrated  in  Section  2.4,  particularly  in  Figure 
2.21.  The  Forgy-Jancey  algorithm  is  known  to  be  capable  of 
finding  clusters  that  bear  no  relation  to  any  natural 
Clustering  [10],  and  must  be  used  with  extreme  care  on  unknown 
data  that  can  not  be  visually  inspected.  Some  of  the 
difficulty  in  locating  the  natural  clusters  in  Figure  2.21  can 
be  attributed  to  the  procedure  used  for  generating  new 
clusters.  Without  this  procedure,  however,  it  is  necessary  to 
specify  the  number  of  clusters  in  advance.  Jancey's  procedure 
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for  monitoring  the  total  summed  deviations  offers  one  method  by 
which  the  number  of  natural  clusters  might  be  determined. 
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3.0  MINIMAL  SPANNING  TREE  ALGORITHM 

3.1  Introduction 

In  this  section  we  describe  the  minimal  spanning  tree 
algorithm  [23  for  clustering  data.  This  graph-theoretic  method 
differs  greatly  from  the  Forgy-Jancey  method  described  in  the 
previous  section.  We  used  an  implementation  of  the  MST 
(minimal  spanning  tree)  algorithm  in  our  investigation  obtained 
in  the  ARTHUR  subroutine  package.  ARTHUR  is  a  general  package 
of  data  manipulation  and  analysis  subroutines  that  was  obtained 
from  the  University  of  Georgia.  The  MST  subroutine  in  the 
package  is  called  TREE  and  requires  use  of  a  distance 
subroutine  DIST,  plus  input  and  output  utility  subroutines  also 
contained  in  the  package.  ARTHUR  will  be  discussed  further  in 
this  section,  along  with  the  description  and  discussion  of 
TREE.  Before  we  examine  how  TREE  clusters  our  data,  we  must 
first  discuss  the  MST  algorithm  and  how  it  can  be  used  for 
clustering. 
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3.2  Minimal  Spanning  Trees 

A  minimal  spanning  tree  is  an  undirected  graph  that 
includes  all  the  given  data  points,  has  no  loops  or  cycles,  and 
has  a  minimal  distance  over  all  the  edges  in  the  tree.  The  MST 
can  use  any  distance  metric  to  describe  how  close  (or  similar), 
the  data  points  are,  but  the  important  point  is  that  the. 
resultant  tree  will  have  the  smallest  total  distance  of  any  of 
the  possible  spanning  trees.  The  choice  of  the  distance  metric 
makes  this  method  general  for  any  type  of  data;  however,  in  our 
study  we  only  considered  euclidean  distance  for  x-y  data.  This 
is  because  we  were  only  doing  x-y  clustering  in  this  phase  of 
our  study.  A  different  form  of  metric  could  be  used  for  some 
other  clustering,  such  as  intensity  clustering.  The  ARTHUR  MST 
algorithm  clusters  data  by  finding  "inconsistent"  edges  in  the 
MST;  the  tree  itself  does  not  contain  clusters,  but  the 
algorithm  points  out  edges  that  when  deleted  would  give  a 
natural  clustering.  This  method  of  forming  clusters  eliminates 
one  of  the  major  problems  encountered  in  using  the  Forgy-Jancey 
algorithm,  namely  that  we  must  tell  it  how  many  clusters  to 
find,  but  the  MST  clustering  algorithm  then  needs  some  form  of 
information  describing  what  an  inconsistant  edge  is.  The 
factors  we  use  to  determine  whether  an  edge  is  inconsistant 
depend  on  the  average  length  of  the  neighboring  edges,  where  we 
can  specify  how  deep  along  the  tree  a  neighboring  edge  is.  We 
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also  need  to  know  how  many  standard  deviations  from  the  average 
the  edge  is  question  is.  The  edge  is  cut  based  on  these 
parameters,  if  the  edge  is  more  than  the  specified  number  of 
standard  deviations  from  the  average. 

This  method  of  forming  clusters  is  general  in  as  much  as 
the  distance  metric  is  general  and  the  clustering  criterion  is 
relative  to  distances  in  the  tree.  This  method  of  choosing- 
clusters  based  on  inconsistant  edges  in  the  minimal  spanning 
tree  will  work  well  on  many  types  of  data,  as  shown  by  Zahn 
[2].  After  the  tree  is  formed,  the  inconsistant  edges  are 
found  and  deleted  and  the  connected  groups  that  are  left  are 
the  clusters.  This  means  that  we  can  form  any  number  of 
clusters,  not  just  a  predetermined  number  or  range  of 
clusters.  Also,  we  can  find  the  tree  just  once,  and  try 
several  different  sets  of  clustering  parameters  on  it,  without 
the  expense  of  finding  the  MST  each  time. 
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3.3  ARTHUR  and  the  MST 

Vie  used  the  TREE  routine  in  ARTHUR  to  cluster  data  using  a 
minimal  spanning  tree.  TREE  can  be  given  three  parameters,  ..the 
DEPTH  of  the  neighbors,  a  FACTOR  used  in  finding 

inconsistencies,  and  a  SPREAD  of  the  number  of  standard 

deviations  allowed  before  an  edge  is  inconsistant.  An  edge  is: 
inconsistant  if  its  length  is  greater  than  FACTOR  times  the 
length  of  the  neighboring  edges,  or  if  it  is  more  than  SPREAD 
standard  deviations  longer  than  the  average  of  the  neighboring 
edges,  where  neighboring  edges  are  all  edges  within  DEPTH  edges  . 
from  the  edge  being  tested  [7].  A  copy  of  the  ARTHUR 
documentation  is  contained  as  an  appendix  in  Fordon  [9].  In 
using  TREE,  the  data  must  be  input  to  ARTHUR  using  one  of  its 

utility  routines.  The  distance  matrix  containing  the  distance 

between  every  pair  of  data  points  must  be  found  and  stored 
using  the  DIST  utility  routine,  and  then  TREE  is  called.  Part 
of  a  test  run  of  ARTHUR  and  TREE  is  shown  in  Figure  3.1* 

showing  how  TREE  outputs  the  MST  as  a  list  of  nodes  and 

neighbors.  The  resulting  clusterings  are  shown  in  Figures  3.2, 

3.3,  and  3.4,  as  lists  of  entities  (they  call  them  patterns) 

belonging  to  the  various  clusters.  These  clusterings  were  all 
done  on  the  same  MST  in  the  same  run  of  ARTHUR;  recalculation 
of  clusters  with  different  parameters  does  not  require 
recalculation  of  the  distance  matrix  or  the  MST.  The  main 
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limitation  is  that  in  examining  N  data  entities,  we  require  on 
2 

the  order  of  N  units  of  storage  for  the  distance  matrix,  and 
of  course  the  amount  of  time  to  compute  all  the  distances. 
This  is  significant  in  that  if  we  double  our  number  of  points, 
we  need  four  times  the  storage  and  computer  time  to  do  the  same- 
processing.  For  example,  a  10  by  10  picture  has  100  points,; 
and  needs  10,000  words  of  storage  for  its  distance  matrix;  if 
we  increase  the  picture  size  to  100  by  100,  we  have  10,000 
points  and  need  100,000,000  words  to  store  the  distance 
matrix.  So  to  enlarge  a  picture  by  a  factor  of  10  requires  a 
4 

10  increase  in  memory. 

A  temporary  limitation  of  TREE  is  that  it  only  prints  a 
description  of  the  MST  in  terms  of  the  nodes  and  edges  as  in 
Figure  3*1 »  and  clustering  cannot  be  displayed  other  than  in  a 
list  of  cluster  members.  The  cluster  membership  is  not  output 
in  an  ARTHUR  format  file  so  we  cannot  display  the  clustering 
using  one  of  ARTHUR'S  plotting  routines. 
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3.4  Actual  Clustering  Using  TREE 

In  order  to  test  the  clustering  ability  of  TREE,  we  used 
data  previously  used  in  the  Forgy-Jancey  testing.  We  took  part 
of  one  of  our  rayleigh  data  files,  namely  cluster  number  2  from 
Figure  2.11,  and  converted  it  to  ARTHUR  data  format.  This' 
picture  consists  of  218  data  entities,  shown  as  1's  in  Figure 

3.5.  We  want  to  isolate  the  obvious  group  of  121  entities  from 
the  surrounding  outlier  entities.  We  ran  ARTHUR  on  this  date, 
found  the  distance  matrix  using  euclidean  distance,  and  called 
TREE.  TREE  found  the  minimal  spanning  tree,  shown  in  Figure 

3.6,  and  clustered  the  data  according  to  several  different  sets 
of  parameters.  The  clustering  formed  by  the  default  parameters 
of  DEPTH=3 ,  FACT0Rs2  and  SPREADsO,  as  in  Figure  3.2,  did 
eliminate  some  of  the  outliers  from  the  desired  cluster,  but 
left  207  entities  in  the  cluster,  which  is  not  very  close  to 
our  desired  121  entities  in  the  cluster.  These  parameters  did 
not  form  many  extra  clusters;  only  five  clusters  were  found, 
seperating  only  11  outlier  entities  from  the  main  cluster. 

By  trying  only  a  few  different  combinations  of  parameters, 
we  made  the  clustering  much  better.  Using  DEPTHsI,  FACT0R=1, 
and  SPREADsI ,  we  found  33  clusters,  shown  in  Figure  3.3.  This 
clustering  put  132  entities  in  the  main  cluster,  eliminating 
almost  all  the  outliers.  It  removed  86  of  a  possible  97 


outlier  entities  from  the  main  cluster.  Using  the  same 
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parameters  as 

before,  but  setting  SPREAD=0.5, 

we 

obtained 

the 

best  results 

as 

shown  in  Figure  3.^. 

We  found 

37 

clusters, 

and 

the  main  cluster 

had  only  two  outlier 

points 

included  in 

it, 

for  a  total  of  123  points  in  that  cluster.  We  were  therefore 
able  to  eliminate  95  of  97  outlier  entities.  This  clustering 
is  drawn  out  in  Figure  3*7,  shown  as  the  remaining  edges  in  the 
tree,  with  all  of  the  clusters  circled  to  identify  them. 

This  clustering  differs  from  possible  clusterings  by  the 
Forgy-Jancey  algorithm  in  several  ways.  First,  it  cannot  split 
up  dense  clusters  because  it  attempts  to  cut  edges  only  where 
they  are  not  dense.  Also,  since  this  is  a  single  pass 
clustering,  the  decisions  are  made  once,  not  iteratively 
updated.  This  means  that  there  are  no  convergence  problems 
with  the  resulting  waste  of  computer  time.  Again,  since  we 
need  not  choose  in  advance  how  many  clusters  to  find,  we 
eliminate  the  problem  of  re-running  the  program  with  different 
numbers  of  clusters  until  an  optimal  number  is  found.  Instead, 
the  program  will  find  the  appropriate  number  of  clusters,  given 
the  parameters  for  cutting  inconsistant  edges. 
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3.5  Conclusions 

The  TREE  routine  can  create  clusters  from  noisy  data,  as 
shown  in  Figure  3.7,  without  problems  such  as  the  splitting  of 
desired  clusters  found  with  the  Forgy-Jancey  clustering. 
However,  TREE  in  its  present  form  is  relatively  expensive  in 
both  time  and  storage  required.  Part  of  this  expense  is  due  to; 
the  overhead  involved  in  using  the  ARTHUR  subroutine  package, 
but  using  ARTHUR  is  worthwile  because  it  gives  us  a  powerful 
and  flexible  system  of  programs,  rather  than  a  stand-alone  MST 
program.  ARTHUR  provides  us  with  its  flexible  input  routine, 
and  has  built-in  routines  to  do  scaling  on  the  data,  and  even 
can  orthagonalize  data  that  is  not  orthagonal.  It  also  has  a 
histogram  plotter  to  check  the  data  distribution,  three 
different  data  plotting  routines,  and  several  classification 
subroutines  besides  TREE,  such  as  a  bayesian  classifier  and  a 
hierarchical  classifier.  A  stand-alone  MST  program  might  be 
more  efficient,  but  would  take  substantial  programming  effort 
and  would  sacrifice  all  those  valuable  subroutines.  The  space 

2 

expense  is  due  to  the  N  size  of  the  distance  matrix,  which 
presently  limits  maximum  size  of  the  data  set  that  can  be 
analyzed.  However,  ARTHUR  can  store  data  in  both  memory  and 
disk  dynamically  by  swapping  parts  of  the  data  set  between  core 
and  disk,  saving  memory  at  the  cost  of  more  disk  I/O 


operations 
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TREE  is  more  flexible  than  Forgy-Jancey  because  it  does 
not  have  to  told  how  many  clusters  to  find.  In  addition,  Zahn 
[23  suggests  a  variety  of  alternative  schemes  for  dividing  a 
minimum  spanning  tree  into  clusters,  each  of  which  might  be 
useful  for  a  different  pattern  of  data  entities.  If  some  of 
these  schemes  were  to  be  incorporated  into  the  MST  subroutine, 
it  would  significantly  enhance  the  usefulness  of  the  minimum 
spanning  tree  as  a  clustering  technique. 


MODELING  OF  SHADOWS  IN  RADAR  CLUTTER 


4.0  PREPROCESSING  TECHNIQUES 


4.1  Introduction 

•In  this  section  we  describe  a  few  preprocessing  techniques 
used  on  the  data  before  applying  them  to  the  clustering; 
algorithms  discussed  in  the  preceding  sections.  The 

preprocessing  techniques  used  are: 

a)  point-to-point  local  variance 

b)  point-to-point  local  gradient 

c)  radial  derivative 

Of  all  these  preprocessing  techniques  the  radial  derivative  is 
the  most  promising  one.  These  techniques  were  applied  to  the 
compressed  digitized  radar  data  provided  to  us  by  RADC  and 
shown  in  Figure  2.19.  In  each  case  we  applied  the  technique 
and  examined  the  output  data  for  clues  for  identification  of 
bright  and  dark  areas.  Results  obtained  for  each  technique  are 
discussed  individually  in  the  following.  For  each  of  these 
techniques  we  used  half  of  the  picture  provided  as  otherwise 
the  amount  of  data  to  be  processed  was  too  large. 
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4.2  Point-to-Point  Local  Variance 

The  idea  behind  this  method  is  that  probably  the  variation 
of  the  intensity  of  the  radar  return  will  be  very  small  inside 
a  shadow  and  it  will  be  large  in  the  peripheries  of  shadows. 
For  each  point  in  the  picture,  variance  was  calculated  by 
considering  nine  neighboring  points  on  the  left-right  and  above 
and  below  this  point  as  shown  below 


y 


XXX 
X  ■  X 
XXX 


We  first  calculated  the  average  intensity  I  for  these 
nine  points.  Then  we  calculated  the  local  variance  by  using 
the  formula 

^  _  2 

Variance  s  L  (I  -I)  /9 
j-1  3 

where  I  indicates  intensity  at  the  j-th  point. 

J 

Then  these  variances  were  plotted  on  an  eight  level  grey 
scale  using  the  following  symbols:  (W,#,A, 1,+,-, . , blank)  in 
descending  order  of  the  value  of  variance  (K  represents  the 
highest  and  blank  represents  the  lowest  values  of  the 

/  / 
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calculated  variances).  This  plot  of  the  variances  is  shown  in 
Figure  4.1. 

An  examination  of  this  figure  reveals  that  the  shadowy 
areas  do  have  small  variances  and  the  edges  have  large 
variances.  However*,  some  areas  with  constant  but  high- 
intensities  also  have  small  variances.  Hence  it  appears  that- 
the  local  variance  alone  is  not  sufficient  for  clear  cut 
idetif ication  of  shadows.  Other  information  like  intensity  and 
edge  information  may  have  to  be  considered  together  with  the 
local  variance.  Further,  with  this  local  variance  method  all 
directional  information  i.e.  the  variation  of  intensity  with 
distance  is  lost. 

For  an  128  x  128  point  data  set  computation  of  local 
variances  needed  33  seconds  of  CPU  time.  The  storage 

2  2 

requirement  was  2  x  128  bytes,  128  bytes  for  the  original 
2 

array  and  128  bytes  for  the  calculated  variances.  We  used 
200  kilobytes  of  storage  and  it  was  sufficient.  Expense  for 
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<1.3  Point-to-Point  Local  Gradient 

The  idea  behind  using  local  gradient  was  that  there  should 
be  higher  variation  of  intensity  on  the  edges  of  shadows  than 
inside  them.  For  each  point  of  the  array  local  gradient  was 
calculated  by  considering  four  neighboring  points  as  shown 


local  gradient  -/<^(x2,yfI(xry)>2+(I(x,y2)-I(x.y1)>: 


where  I  represents  the  intensity  at  the  point  (x,y).  The 

(x»y) 

calculated  local  gradients  were  plotted  over  an  eight  level 
grey  scale  using  the  same  symbols  as  in  the  case  of  local 
variance.  This  plot  is  shown  in  Figure  4.2. 

An  examination  of  Figure  4.2  shows  that  in  the  case  of 
local  gradient  also  shadows  and  high  constant  intensity  areas 
have  low  gradients.  The  edges,  however,  showed  larger 


gradient.  Hence  the  local  gradient  above  also  can  not  be  used 
to  identify  shadows  conclusively.  Other  information  like  the 
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intensity,  needs  to  be  considered  together  with  gradient. 

For  a  128  x  128  array  of  data  points  it  took  25  seconds  of 
CPU  time  to  compute  gradients.  Storage  requirement  was  also  2 
2 

x  128  bytes.  An  allocation  of  200  kilobytes  was 
sufficient.  Expense  for  one  run  was  $1.91. 
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4.4  Radial  Derivative  Method 

In  this  method  derivatives  were  calculated  along  radial 
lines  along  those  of  the  radar  beam  starting  from  the  location 
of  the  radar  as  shown  by  the  X  in  Figures  4.3  and  4.4.  The 
idea  behind  using  radial  derivatives  is  that  as  the  radar  beam 
travels  in  radial  directions,  the  variation  of  intensity  in  the 
radial  direction  will  be  a  better  measure  of  the  reflective 
nature  of  the  objects  illuminated  by  the  radar.  Also  any 
object  intercepting  the  radar  beam  will  produce  a  shadow  in  the 
radial  direction. 

For  each  point  of  the  data  array  radial  derivative  was 
calculated  by  considering  four  neighboring  points  as  shown 


<I(*.y2rI(x.y1))^l',f(x2.y)“I(* 

Radial  Derivative  ■  - - - - = - 

/z<x2+y2) 


for  the  points  in  the  first  quadrant.  This  formula  has  to  be 
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modified  for  other  quadrants. 
For  example  for  Che  points  in  the 

fourth  quadrant  as  shown,  the 
radial  derivative  is  given  by 


Radial  Derivative  *  ^(x,y^) 


(x,y2)^y^+(I(x2,y)"I(x1,y)^x^ 

V2  (xV) 


These  radial  derivatives  were  plotted  on  a  seven  level 
grey  scale  using  the  following  symbols. 

#  <  -30 
-30  =<  Q  <  -15 
-15  =<  -  <  -  5 
-  5  =<  Blank  <  5 
5  =<  +  <  15 
15  =<  A  <  30 
30  =<  W 

the  plots  are  shown  in  Figures  4.5  and  4.6. 

An  examination  of  Fig.  4.5  reveals  that: 

1.  All  dark  regions  have  low  radial  derivative 
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2.  Range  marker  showed  up. 

3.  Shadows  were  marked  by  the  following  characteristics:  in 
front  of  the  shadows  we  observed  regions  of  high  radial 
derivatives  (cause  of  shadow)  followed  by  regions  of  low  radial 
derivatives  (edges  of  shadow)  further  followed  by  regions  of 
lower  radial  derivatives  (shadow).  It  was  also  observed  that 
shadows  were  followed  by  regions  of  higher  (increasing  radial 
derivative  (elevated  region  beyond  the  shadow).  We  also 
observed  regions  of  negative  radial  derivative  before  shadows 
and  regions  of  positive  radial  derivatives  after  shadows. 

4.  For  uniform,  but  high  intensity  regions  (non-shadows)  we 
observed  regions  of  positive  radial  derivatives  before  them  and 
regions  of  negative  radial  derivative  after  them.  Of  course, 
these  uniform  intensity  areas  themselves  have  low  radial 
derivative  like  shadows.  Thus  the  signs  of  radial  derivative 
before  and  after  a  region  of  low  radial  derivative  may  be 
useful  in  distinguishing  shadows  from  other  low  derivative 
regions. 

5.  However,  we  found  some  exceptions  to  these  observations 
also.  Radial  derivative  seems  to  be  the  most  promising  of  all 
the  preprocessing  techniques  test  by  us  for  detection  and 
identifying  shadows. 

However,  in  order  to  put  confidence  into  the  use  of  radial 
derivative  for  detection  and  identification  of  shadows,  more 
testing  with  data  with  a  corresponding  topographical  map  is 
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necessary.  VJe  did  not  have  the  topographical  map  to  verify  our 
conclusions.  Also,  the  radial  derivative  makes  use  of  radial 
changes  and  it  does  not  take  into  account  of  side-to-side 
changes.  Using  both  radial  and  local  derivatives  may  improve 
the  decision  making. 

Originally  we  plotted  the  radial  derivative  information  or 
a  seven  point  grey  scale.  However  this  plot  looked  too 
cluttered  with  symbols.  Hence  to  have  a  clearer  picture  we 
replotted  the  radial  derivatives  on  a  five  level  grey  scale. 
The  derivative  values  represented  by  +  ,  blank,  and  -  were 
replaced  by  blanks.  This  plot  is  shown  for  the  upper  quadrant 
in  Figure  4.7.  The  regions  of  low  intensity  depicted  in  Figure 
2.20  are  shown  dotted  in  Figure  4.7.  Using  the  plot  we 
examined  the  regions  of  high  but  relatively  constant  intensity 
preceeding  shadows  to  test  them  as  precursors  of  shadows.  We 
noticed  that  for  these  regions  of  high  but  constant  intensity 
before  shadows  there  were  regions  of  positive  derivative  before 
them  and  regions  of  negative  derivatives  after  them.  This 
indicated  that  these  regions  are  some  sort  of  precursors  of 
shadows.  (These  areas  were  selected  by  hand.  This  selection 
may  be  done  by  a  computer  program.  However  these  programs  may 
be  quite  expensive  and  it  is  easier  to  select  them  by  the 
sight.)  Figure  4.8  shows  the  same  information  as  Figure  4.7, 
but  with  the  precursors  indicated  by  numbers  0-9  and  letters 
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B-U. 

Then  we  tried  to  distinguish  these  groups  of  bright 
regions  to  see  if  these  belong  to  more  than  one.  type  of 
regions.  Our  investigation  did  not  bring  out  any  significant 
differences  among  the  various  groups. 

Next  we  attempted  to  correlate  these  precursor  regions  to 
the  shadows.  We  did  observe  that  there  were  shadows  following 
many  of  these  precursors.  This  gave  confidence  to  our 
conclusions. 

Finally,  we  used  the  spanning  tree  algorithm  tc  cluster 
the  dark  areas  (possible  shadows).  We  did  not  obtain  any 
accurate  division  of  these  areas  into  individual  shadows. 

For  a  128  x  128  point  picture  it  took  26  seconds  of  CPU 
time  to  calculate  the  radial  derivatives.  Storage  requirement 
2 

was  2  x  128  bytes.  We  were  alloted  200  kilobytes  and  it  was 
sufficient.  Expense  was  $1.95  per  run. 
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4.5  Conclusions 

Out  of  the  three  techniques  of  preprocessing  used,  the 
radial  derivative  technique  seems  to  be  the  most  promising  for 
detection  and  identification  of  shadows.  However,  combinations 
of  these  techniques  may  prove  to  be  more  useful  under  some 
conditions.  It  is  necessary  to  have  a  topographical  map> 
corresponding  to  the  data  used  to  verify  our  conclusions. 


! 
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t 

4.0  Note 

i 

We,  at  Clarkson,  wrote,  debugged  and  tested  the  programs 
used  for  the  grey-scale  plots,  calculation  of  local  variances, 
local  gradient  and  radial  derivatives.  These  programs  are 
specifically  written  for  these  types  of  data  (type  of  numbers, 
data  x-y,  rather  than  radial  format).  These  programs  are  now 
available  on  disk  and  also  on  tapes  for  further  use  on  this 
type  of  project. 
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5.0  CONCLUSIONS 


This  report  has  presented  the  results  of  a  comparison  of 
nonparametric  clustering  techniques  and  their  use  in  examining 
radar  clutter  from  shadowed  terrain.  Successful  techniques  for 
extracting  features  from  clutter  will  aid  in  understanding  the 
underlying  structure  of  the  clutter  and  in  turn  wiljl  lead  to 
the  development  of  useful  clutter  models.  The  Forgy-Jancey  and 
minimum  spanning  tree  algorithms  were  tested  on  both  simulated 
and  actual  radar  data.  Some  preprocessing  techniques  were  also 
tested  on  the  actual  radar  data  in  connection  with  the  minimun 
spanning  tree  algorithm.  The  Forgy-Jancey  algorithm  was  found 
to  be  a  very  flexible  tool,  handling  large  multidimensional 
data  sets  with  minimal  cost  in  CPU  time  and  storage. 
Unfortunately  the  performance  of  this  algorithm  was  found  to  be 
inconsistant  in  terms  of  finding  the  "natural"  number  of 
clusters.  In  general,  repeated  applications  of  this  algorithm 
are  necessary  in  order  to  find  the  optimum  number  of  clusters 

and  parameter  values.  A  more  detailed  discussion  of  the 

conclusions  pertaining  to  the  Forgy-Jancey  algorithm  can  be 
found  in  Section  2.5. 

The  minimum  spanning  tree  algorithm  was  found  to  be  more 

robust  in  its  ability  to  detect  a  natural  clustering;  it  was 
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never  observed  to  split  a  natural  cluster  as  was  the 

Forgy-Jancey  algorithm.  The  ARTHUR  implementation  of  the 
minimum  spanning  tree  provides  flexibility  in  clustering,  since 
once  the  tree  is  found,  the  data  can  be  clustered  repeatedly 
without  recomputing  the  tree.  Finding  the  tree,  nov/ever, 
turned  out  to  be  costly  for  large  data  sets,  since  the  ARTHUR 
implementation  requires  calculation  and  storage  of  the  distance 
between  every  pair  of  data  entities.  Alternative  procedures 
for  computing  the  tree  are  available,  however,  as  are 

alternative  procedures  for  clustering  the  data  from  the  tree. 

A  detailed  discussion  is  presented  in  Section  3.5. 

A  discussion  of  the  preprocessing  techniques  examined  is 
presented  in  Section  4,  with  the  conclusions  in  Section  4.5. 
The  preprocessing  was  motivated  by  the  need  to  reduce  the 

amount  of  data  provided  to  the  minimum  spanning  tree  algorithm, 
but  was  found  to  be  a  useful  technique  in  its  own  right  for 

extracting  features  from  clutter.  Computation  of  the  radial 
derivative  was  the  most  revealing  of  the  techniques  studied  for 
locating  and  defining  the  edges  of  shadows. 

Clustering  techniques  appear  to  hold  great  promise  for 
examination  of  the  structure  of  clutter  returns,  provided  that 
the  techniques  are  applied  with  some  care  and  common  sense/  A 
primary  result  of  the  work  to  date  has  been  installation  of  the 

*  f 

relevant  programs  on  the  Clarkson  College  computer,  and  the 
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If 

» 


familiarization  of  the  researchers  with  their  use.  The  next 

i 

step  would  be  a  detailed  application  of  these  programs  to 
actual  radar  data,  including  an  examination  of  a  wider  variety 
of  clutter  properties. 
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APPENDIX:  PROGRAM  LISTINGS 

This  appendix  contains  the  listings  of  the  programs 
developed  and  used  in  this  project.  Figure  A.1  is  the  comments 
section  from  subroutine  EXEC,  and  describes  the  input  format 
for  Forgy-Jancey  data.  Figure  A. 2  is  the  input  specification 
for  the  outlier  Forgy-Jancey  subroutine  EXEC3.  Figure  A. 3  is 
the  subroutine  EXEC  to  go  along  with  the  comments  in  Figure 
A.1.  Figure  A. 4  is  subroutine  EXEC2  which  allows  outputting  of 
the  sorted  data.  Figure  A. 5  is  subroutine  EXEC3,  the  outlier 
calling  subroutine.  Figure  A. 6  (Parts  1-3)  is  subroutine 
KMEAN,  the  original  Forgy-Jancey  clustering  routine.  Figure 
A. 7  (Parts  1-3)  is  outlier  subroutine  KMEAN2.  Figure  A. 8  is 
the  RESULT  subroutine,  which  outputs  the  results  of  KMEAN. 
Figure  A. 9  is  subroutine  RSULT2,  which  allows  output  of  the 
sorted  data  array.  Figure  A. 10  shows  the  driver  routine  MAIN, 
function  DIST,  and  subroutines  USER  and  USROUT  for  a  sample 
Forgy-Jancey  clustering  run.  Figure  A. 11  shows  MAIN  and  DIST 
for  a  sample  outlier  run,  and  also  another  version  of  DIST  used 
in  intensity  clustering.  Figure  A. 12  is  subroutine  RAY,  our 
rayleigh  random  number  generator.  Figure  A. 13  is  the  program 
that  created  the  rayleigh  data  set  with  sigma  s  10,000.  Figure 
A. 14  is  the  program  that  created  the  rayleigh  data  set  with 
sigma  =  50.  Figure  A. 15  is  the  plotting  program  that  produced 
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all  of  our  100  x  100  and  128  x  128  point  cluster  plots.  Figure 
A. 16  is  the  grey  scale  plotting  program  used  to  generate  all 
grey  scale  plots  without  overprinting.  These  programs  were  all 
either  written  at  Clarkson,  or  installed  and  modified  at 
Clarkson  to  suit  the  purposes  of  the  project.  For  a 
description  of  ARTHUR,  see  Section  3»  and  the  ARTHUR  program' 
documentation  [7]. 
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Figure  2.1:  The  Forgy-Jancey  Algorithm 
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Figure  2.3:  Test  Scene  for  Clustering 
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Figure  2.4:  Intensity  Clustering  of  Test  Scene 
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Figure  2.5:  x-y  Clustering  of  Test  Scene 
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Figure  2.6:  Intensity  Clustering  of  Simulated  Shadovts 


MODELING  OF  SHADOUS  IN  RADAR  CLUTTER 
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Figure  2.8:  x-y  Clustering  of  Simulated  Shadows 
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Figure  2.10:  Low  Background  Level  Run  to  Convergence 
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Figure  2.11: 


x-y  Clustering,,  Low  Background  Level;  NC=2 
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Figure  2. 14:  x-y  Clustering,  Lov;  Background  Level;  NCs5 
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Figure  2.15:  x-y  Clustering,  Low  Background  Level;  NC=6 
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Figure  2.16:  x-y  Clustering,  Low  Background  Level;  NC=7 
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Figure  2.18:  Test  Scene  With  Outliers 
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Figure  3.1:  TREE  Example  Run 
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Figure  3.3:  Better  Clustering 
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Figure  3. Best  Clustering  Found 
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Figure  3.5:  Data  Set  for  Tree  Example 
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Figure  3.7:  Best  Clustering  Found 
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Figure  4.3:  Location  of  Radar  in  Upper  Quadrant 
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Figure  H.5:  Radial  Derivative,  Upper  Quadrant 


MODELING  OF  SHADOWS  IN  RADAR  CLUTTER 


«•  «***4aa« *%• •*•*••« a  **%•*«•••••«»*• 

tUA*  ••  **••**••  •  •• 

««  '  M«U 

••A*  *  •••  AA#  •  ••  till  •*•!»»•  tl*  <*■•«»" 

-•••iaa-*  44Aa-a-  —  a- -  •  *  a  ♦ 

•#««•  •  —  -••«  itti  —tit  •••• 

-a*4  a  4*.  ••  •  •  •  **•*  ••  *»**  •«••• 


tlU||t|t|M»»(k»|(*r|«t4U 

* »  I  • a • t  •  •  «  I  • A* •  •  •  •  •*•  •  ••••» 

•  •  •  lit •  »•  ♦  •• 

•  I  - Aa**4-a.4*a A*--- a-  •••  »• 

•  *  •  a  a  1441  IttAf 

—  .  •***•*♦•  •*»**•»* 


..•«*  — -• 

•  Mttt  ••  «  •  •  I 

•-  M  —  I44IM 
•••  •*•••  till 

--  »  -  444 


*•  •••!•  4t«4>  ••  #•#•  HIM  ««••«  «l444i»»t»»**tt»t*  44  ........4*. 

•  ••  I  •••4  44  -  -  44  4  #4*. a  A  4  4  a*  4  4  4  •4ll4i*4ll  •(iiti  t  itl.fl  Mt  M4«iMt«l  >44**  —  •  •  A*  Ii44  t4  I  »•  »  4  M  4»  I  4*—*»  ...  4* 

•«  44  44—  44  44  *  *  •  •  a*  *  *  44  4  . . .  .  )l  l»,  «.4I  «4  4  i.4>>  /•••  At*— *4  •  4  4  4  4  4  •  4  4  4  •  4  •  *  --  -t  -  *  % 

«M>  4444«4'4f  441  UIAA4  -  >4444  f  .'«l  444444.^4  1.4.  •  44  .44  *4  .  4  4  4  4  4  A  4  4  4  -•>•  •  4  *  #44  •  4  I— —  ••  4  4  4  ••••  «#  4 

|4C  A4A4I4  •« 44  4 •  44  44444  •.«•  4 *4  ••••••,  |4  4  4  4« 4 •«•>•%• 44  4  *44  4  »•••»•  4»4>t 1 ■• it  ••••4<444i4*-*  4*44  444  -----  ,/ 

«A<4«444t4  f  •••—-•  4  44*4**44  -•«*#  *  ••  J*  •••*  —  44  4  4414  44  *•  4  4  a  *  aa  Att  a  •  •- 

A4^>i4ll*4  •  >  ••••••••  4**444*  4*4  *  A  4  at*t  >»«4*4  44  -a  4414*4(4  444*4  *  •*•4(44*4  .  #  -  ••• 

»«  144*.  4  •>••  -  *•**••*•»«  4**4  >■*••>  —  4**4  *4*  4  4  44444  44*44  — •  ■•J-I4»I4*4*444  ••  4  •• 

*44  44  —  4  •  444**4*  44  444  -  •••!**»  ••  4*44**  44  >•>•*#•>  4444  *4*4  4  ••••(•I  Hi*4(  *4*4  -•*'•« 

4*4444  44444444  •  *4*444  4444  44.44a. .*4  — I.-  4*4  •-  (AAAIAA*  •--.#•*  *•••*••*•  >4..  44444«*4  ••44114 

•  4*44  4  4  *•  ill.  4  4***4  *444**4  4i«l\|lll  •••(*•  4*4  •**••  4**4  —*li  44  •  -  4  4  4  A*  *-*•  *4  •  *  •  ’**■  *•>•**•• 

•*«•  14444  4  14  1444  (AAA  4  44  4  4  44  4  44  •**••  **  4  •(##»*•  44  ••’!•  44  4  AAA  A  A*  a  A  4  -  ••••«*  *.('  A4a 

■  M>  4  4  *  \  4 ->  ••••  44  4  4.  '4*4  444  4  44  »•••  A*  A  4  --til  •••  4*4  •  *-  4*4  «  44444****  •  ••4*‘4i**l4t  •• 

•  ••  |»  44  <;#•  •-!•'•  4***4  4  4**444  —4*  tta*  .*44  «•*  •«  -  ****4*44  4444  -•-  •*••• 

44  *444-  ••••  444  4***4  -  ••  4  *  4  -fa'  *»*•---  •••!  ****** *  -  44  •••  •**•  4  4-444 

444  —  4  44  I  >  4*44  -  -  ******  •  -('*-*•  4  *  •  4  *4*4--  I  A*  *  4  44*  *4  4**A  ‘  * 

•  4  ----  •**•-••*-  —  --  ****44*4  —••»•••  -  44*44  *4  t4-4((*tl  4***4  V*  lt»>  44  *4-  4  444  •**••*•  •• -a  *•••  .  **  •“ 

•  *4  -4-  4  *44  4  J44.l--i««  4  •  at  .1  k  *  •  --*•!«  — 4-  **4*4  -•*-  *(4**4**a*4-  l*4i«44  *4*  --  ***•».*•  -  •••.title* 

•  *l»4>-  4  A*  4 . 4  -  a  t-  1  •  4  t- \  1  •  •  a  ••  #ti  —  • 4lt*4  >>ii  *••»•**••  «4*-4t-ti*#  4*4*4  t«Aktt»l4»(  -t.4i«  — 

•  ••  -•  A—  4«IM*Jtt  4  —  *4  —  44  I  •*  a  4*  4.  *  44—  4 1  *  4*  4  -  •  4 »-  4  4  4  4*  *  •  --444*  *  »V-  4  *  *•  4  *4 «  *  4  4*44  •••'** 

II  •  >|.  4A*4*  — 44  — *444*  «4  ••4tl4><tt«it-  44 4 4 *••*•  —  t «4(—  4*44  *4  4  4  4  4/.*  --  4  •  -  • 

■  A-  -4-  4A*  HI  4  «||l4.i4tlll(llll^-«  —4*4  4*4-(4tt4t'a«i  *444t**4*4  *****  -4-  A* 

HA— ....  *4  4  4*2  «t#4  t«tttil44«A.lf** — —  —•<  «  *  •»MlltH(\4»W»M*Vt  A  At  a  »  ♦  --  -  — 

•t-.*.i«4i.  —  —4*4*  — — 1<  4.  — — — -—  ••••»•>•— ill  t|  -  444  ••(••ill  *****  4***4  —  4**44*11  •••«•• 

—  — A*.  •(-  —Ilf  -—  4444  -  —•*••--  4**44  -*44a  •••  4  ••  •*••*»•••  — 

*«  —lit  4*  -*»•*-  —  4**44  -  *  *  -#••*'-  •4***4444444  -  - 

|A*  ••!«•  44444  — —  4*4444  — — -  •*  *4**444  *tltt  -  *4  *4  444  - 

•  4*4  -4*-  A*«44  -  — —  -  444  — -  4  4  4**4**  -tt  *  t*4  44  -  Jt  A-  — —  *AAAI  --—  *4 

*.  «»  v|.  44*4444  —44*-  44*4  —  4  A***  f  *4*  4.4  at  AAA  *  —  4  -  4*44*44  —  — 

■I*  44*4  444  —4—  444444  -  A*  4  %A  I  A*  ***  4*4  41  *  *  »  -  •  —  —  -—  —  —  .  4A#4»4*4  -  - 

— |A—  4  4***  4  444444  4  At  4  **  *A  4  4  4  4  A  •  4  4***  4  4  “AillAI"  llatlAII  -  444 

4 _ -4A4  —  4«*4  44  •— .  4444  44*4**44  IAA44I4  4  «  — *4  44  4  4  4i  i  4  4  4  •  «' 14  *4-  •*»•  44  •«**• 

•  •  -  *4*4  4**4  —  444  --  4*44  44  •  *♦  -» f  44  4  44  *  j-4  4  A  4 .  i  4  4  _•  -  — -  •»•►» 

--•***  4  *A*  --  —  44  4-*-  —  —  4***4  -44  t-iiAA-tl  •  ••#-  —  -  --  •'  i 

-  4  444*  -  —  —  »IAAt«—  4*444*444  -  -  J  t—  —  4a  *(»•»• 

—  -M*ti-  44**4  #4  —  *444  -<•  4  4  -*44444£(  ••IMiAiHI  - — -4.l'4il>4 

.|»l|  444**44  4444  -  4  *A*  A  —  4*4  Ittltttf*  *4Mli|t*444  —  —  ••italic 

..|t«.  44**44  4*4  •»»**  (tAAAl  -»•••*-  4AA*t»*  *4444444  —  —  -  —  -  4  4  *» 

•|b—  4*4*44  444  -—*****4  4  *****  ♦  -9*-  44*4*4  *44  —  44  -•*- 

_ *_  4*44*44  4**44  -A—  —  4At«(AA4  -*f-  4*44*4  444  4  4  —  -  444 

4  -«*-  **4*4*4  4*4  44* « 4—  a  *  a  al  —  ’  (>■  444*4**4  -til-  **444*  *  -  4  444  -  444  —  A- 

4  4  ••|«|—  4*44*  4  4*44*4#  AA-44  4  4444*4-  4*444*4  4**4—  4*44*44  —  -  4  ••*•••  4  •**•••--  — 


4 4**44*444*4-.. ttillA 
4  44  4  **  ** A •**  —  - 

—  4*1 -  •*«••***** 

•Atll—  4*44*44 

-•*•*—  4**4  **- 

-  - 4  44  4  J-  44 

—  -Aft 4  #4--  --  4*4  —  - 

-A444AA-  -4-4**4 

— «***••-  —  *—  4*4 

— -lift «.44  —  4 

—..*•144-44 
---#••••  B 
4  -  —  (144444a- 
4  — -  ilHIU* 

•••Mini*  - 

—-**••*—— 

4  —  a|#44i-«#— 

4  — .|*t|.ai4| 

•  -  — •••aiiAa- 


-WM- 

4  «444i- 

•A*  Atifl- 
4 A#*  *444  1- 
4*44*  illtt- 
•  ##A*  -•••(- 
—*•4*4*4  •t«l»' 
•4(14  **44  -«»i - 

•4—*44  •  »4*B” 


4  4*44  -ill-  4*44  4  4 4A*4»44 44*4-* .4 it\ JA  4***4' 

•44  — •**»  4  44  4  **  **  *  •  A  A -  —  ♦*•< 

-  -  —••«.  4***4  —  44444***44**444  44  - 

— -  -tiff-  44**4  —  4A«— -  4*4  44***44  444*. 

-Alt —  4***4  -At**--  4*44*44  44*i 

444  •  —  I-  A*  •  A  -  till—  4**4  ♦*•  4*1 

A4AAA4  44*44  4*4*4  - - A  44  4  I-  44  >**■  *1 

*A*4  ********  4*44*  —  -A»t4t4—  —  4*4—  4  41 

44*444*4**4  *-  «  4*4#  *  -«*#4*A-  -*-4**4  •  « 

I  44  *****  -*#**-  4*4*4  —  ««•***-  -A-4A4  < 

AAA  •— - iltllli-  4*4*44  —*#111-14  —  4  — 

ltaWt44444  4444«  4**4. 4  — *44  44- 44  -WM* 

**4*«4A*44l 4 A • A  — —  4*44*4  *•*•»•••■  *  «****• 

|44IW4t— -- #4<  A-  4  4*4*4  4  -  «—  444444a-  AAA  4*444- 

-  —*4 1-  4*44*  4  —  444*1  *A-  4*4 A  <***)- 

A*  -.4*  4*4*4  —-*#*•*»•  -  4*4#*  A4  444- 

**4  *44  -tit-  141(11  ——  —  I  A*  •*————  44*4*  -444*. 

•  44  444  -4*4*  *4*44*4  ♦  — *#t*A-4^-  — 1-4*4*4  -tali 

•  •  A*  4>-  4***  4  —  •*•#-•*«*  -4(14  A*4  4  #i  - 

f  A  —  -4*4-  •«••«  •  -  —  itt.Uta-  -t.tii  •  -*•#- 

4 A-  •  .»•*-  4*4*4  4  -  •  -*«-  - A  A  J  A* •  -•-  *  *4-44* 

A .  -  -  —  —  A  4  *  4***4  4  ♦*  ♦  --  -A  #•  A*  —  A  t  #  A#  4  A  A#  *  *  4*  |ti> 

-  44  -*#t-  *4*444  4  4  —  4  •*«  *#  At A*4«*4# *4 4*A 4*  #4— 

4444  4AA#«  -*--  •  4#  A  *  4  44  —  44AAA-A-  44  •**«  4  4  A « 4- .  44- 1 

••****44  44*4*4————  4*4*4  —  44  -  — -  4***4- 4  —  444AAA4  — -• 

||4lt«AA4  - 1 4J-  444*1*4  -  *44  444  -  ••***«  -  44  444AAA4  •* 

| *4  444*44  -<|«A44#<  ••••*•  -  «*%*«*«**  444**4  44  4444444*4  * 

•  *  *A  *•  t*4  -Mi-ittr*  4  *•#  A4  -  4**4*  A  44  4  4444  4444*««#* 

•  444**4  — —  4***4  •%•!•  --  *44**4  4  44444A4A* 

-<4C-4«»kt  At  M  -.444-  •** « *  *444*4  --  *44  44444*44* 

**•••*4**4  — %  *t  —  4  *  A  A  4AA4A4A4  -  -  44444**** 

•4444  4A»«4*A*A4  -< I-  •*#*•  *A%#t44*4  -  44444**4 

4A44  *4*4  44*  -•-  4 A#  t •  •**■•***•  —  4  4 AAA 

4 A* •  *44  4***4  4*4  *4  4 •  —  4*4*4#  —  -  444A 

4  AAA  4  4*4  44  4444*  •*•--  -  *4*4  -1*4-  -  —  *4* 

•  4  4*|444*4*4  -  4  #4  —  • • A—  A A  4  ••»•••  *4  -A#*--  44 

-  4  AA  44*4  4  •  A  ii.  •’44*44  lIMI—  #*#*-#AtB- 

—  44444A.4*  —  t «  44-  --4*1*444-  4AA4  «•••*- 

4  * A*  *A4  444*  — — —  —  —4 (•**•*—  4 AAA  4***4 

•  44  4  AAt •*•  A»44»A*  --  -4*4***-  — 4AA4-444  4- 

•  14  4AA*#4  44  #4444  #  A  --Ail  •*-  44  —  AAI-il4|— 

44A4  4  #  4*  A  •  A  -  44*4  -it  AA*-t*4»— 

4  4  4  4  4  --  •*•«  4***4  —  w#t  4  A  A  -•**!- 

••  **##  AlAAi*  -***•—--  44AAA4  —  4t-4A4  tail* 

III*.  — — —  4AI4A444  AA*  -  4  4lt--(«i -  4*441  -4 A-4  *4-*4t- 

I4fl*  •■*#>  4  A  A  At  tl  #4  4*4  •*  4—4  •#•- -  •  *•  #  4  -I-  *  -•#•#- 

a-  —4-  •  444*44**4 --  •*##•••  -  •  *44*  •*•  -*•«•*- 

•  AA#  4  A  4  •  44  —44  4444*4  4  4444  |-  -  4*4*  -  — ■  -«Af*BB  — 

— —  •  AA*  *4  A  4  4AAAAAAA4  *4  . 4  -  -ii—  AAA  4—  4*  -*-  -*%*•#- 

f  444*  444*4  *  *****  ***  -»-•  -•  -A###  AA*  —  *•  -A-  *  -*•#•- 

•  f||aa  «  **•*•  *A«  -i*  ---  —  at-44*  —  •**  --  4  *•••»• 

•  4*4*4  •#*!-«•  —  -i-  4*  -t  AA**---  4  -(••*- 

A  A**4  #AA*  AAA*  -  4*  -  *  —  4  14-4*14  --•  *  «••*- 

•  •#  ft  *****  44*  •  444  — -  44-  *  A A*  -a-  *  -#**• 


•  -  •  —44—  -it J  AA*  ---  4  44—441 

44  4  --  -A  *4  Ai  —  A  4 1  A*  4  A  A#  A  4  *A 

44  —4  AAA  •#  A* AA4 A At# A4 4AA •• It— 

*4  —  44AAA- A—  *4  4AAA44 A44-.A*— I—  4  — 

♦  4  -  —  4AAA4-—  4  —  *44  4*4  4  -  — 

•••  444  -  4 • AA A  4  -  44  444AAA4  -•**•#- 

AA%AA A A*  *  444AA4  44  444**44*4  4#44|- 

•  44*A A *4  4  ••**  444  4*4  A#A  It  44- 

--  *••**•  4  4*4*4*44*— I*  4A— — 

•  44444  ••  •••  44 444  A  44 A* - >4. «-• 

•  AAt  A*  44  - -  44  44  4  At  A*  A  - I4lf* 

**44444*4  -  44444 AAA*  --  “ABA  — 

•  **■•**#•  —  4  4 ***** - #»#I* 


•  4  -#•*-- 
4444  — **4  I  — 

4AA4  4# #•• — 

•  A* 4  4*  44 A 
— ***•-•*• #- 
44  —  A Al  - 1 1 •  I— 

444*4  -alt*#*  —a  #tl  — 


4AAA4444  *««•  ******  ••———-  —  •  •  •  a  •  A*  AA  •  •  4,  -- 

4444*4  -*••«)  **A4*A*4  -----  44  4  A  AAt  a  A  *  - 

#4  4  A *4  (•••*4**44**444  —  A*  4  4  4A  A \ 4 

444**444  -if#*  AA*44  444A*4-il  1  4  • 

44AAA44*  —  ta  ti  4* A4  AAA*  -.1# 

4AMAAA*  —AAA  A  A*—  AAA**  4  <  «- 
•A»4*A4  »S#-4  *•-*  *****  4— — — — 

4  •**###*•  -tt-4*  •“**#*• *  *  •- 
•  44  At— WA*  .4-  -44  A  ***** 

♦  4A##«*A«  - C444*  ••*  4 

44  ****■%  A  *  4  '141*1*1- 

•**#•#«••  —it •••••;—— 

44*1***444  -|U#f  •••  — 

4«*t(**4  4*  —*»#•»*.*  - 

44*t#t#A  44  -ta  — - 

44A At  t#A  *  4*  't*»»l«/»4*' 

******  I  44  - 1 »  4  •■«•«•  — 

•  4  Afett#A4  4  Hall"  -- 

—  •**•#•*•  •  -..*•»(.«•*• 

4  *# •#  A4  444  -•11111.111' 


»  ••  *•#•**  -  444  ••4#«a.ati*  — 

■’  -  4  4*44  44  4  *•  —  —  tt  .(•*••* 

t-  —  4  44#  Vt*  A  -  ••#  — --ti-A.t 

fct-i-  4  —  •**•••*  4  •  •#•  —  —  —  a  •  #  1' 

—  —  AA-  4  4*4  lit*  4*  4  --•••* 

-••••A-  4  At  #••**•  •  44  -tt.. 

A4  t # #  9k—  44«lalM*  •  —  t  a  *t 

*•*  me  ••#**••«•  --  -  i  *  •  a 

At  a  *-  ll»i—  •♦*♦».*»*  -  —  -'  •» 

*••**- .4.  W--  •*■#•* 

A4AAA  - •(*(-  44 *•••*•-—-- •  .  a 

*****  --  -ABA—  —  1144***4  •«**# 

•  At  A  A*  —-••#{-  —  -  #  #  At*  » •  4  4  4  4t 

•  •**«*•  -i;-  4  —  -  4  *#t • A#  —  t -4 

•  *****  •*•”  — 1(—  •*••■••  *4 

44**4  44*4  -tt •-  **•*•••  - 

4**4  • AAt  •##••  —  *  4*  *A  t* 

•44  44  •**•-••••-  — -  4*44*4 


4441**444  '#«•*  4  --f- 
•4444**4*44  -*•»#-  -at. A 


4  4  AAA  4  4.— 

•  4**4  *••- 


tt  AA*  -••««'-(•(•  •*••*  -ti-4  44-t4t-  ••***«**  ••«• 

*  #4  «A4  *| t--t  4#t- -  **••*  4  -alai-  ********  — •*»'•*.' 

>tAt#t*A*--  -••#>--  -  -  At#  A  -t-  -tttti-  4  4AAAAAA  *  — •*  .4--  - 

4  A  •!••*•  4  1 4  •  4  I-  -  -•-  •  A  •  A  -  ----  -II  <4 1  J— —  4  4  4  A  AA****  --a*  a  4  «  -  — — - 

4  *******  •  ••  .at--  -  -ti-  AAA  4—  44  -I-  -l(44(-  444AAAAAA*  A  ••  --##•- 

•  4***4  •  *  A  --  -A  tt#  A  A  *  —  4*  -  #-  4  -  4  •••-  44  4AAAA*  *  "*»•*-  •• 

1*4*4  •*•  -4-  —  -4-4*4  —  •*•  --  4  -•••*-  4  •  AAA  *  *  •  •  *  -  -  a.  a  *  *••  He 

4  At  *4  *  I  A*  — —  •  —  -i-  44  -t  AA**---  4  -*#•*-  4  •  A#  A  *  444*  -  *  J—  •  t  *t 

4**4  •  AA*  -  44  •  4  — •  4  («-•**•  --•  4  at  *t-  •****•  *4*  -  -  -  - 

...  ft  4  AA*  4  •••  -  ••*  — -  tt-  *  A  A*  -•-  4  --*■•-  4AA-A4  »••  ----»••*•  -tt* 

|a.  •  — --  4  *  A  A  A  A  *•  *  -  444  —  -  — —  4*  4  -I-  4  •»(>•  *****  •*•  ---tat-  -t  -  - 

|1>M.  ••  —  •  44  --  ••  •••••**•  —  **»  —  —  tt  4*4  ---  4  -ttttti-  4  At#  A  4*4  - ----- 

•4aAa  aaMaa-aa.a— ••  —  ••••  -•  44*  MU  —  *•••  —  — -  444  ---  4  -tttttt#  ******  -  •••  -  —  -*■*•>  -  — 

JT'Tt— —  a|||4— — .  •••  444  —  44  <  A  ■  •  <  -.-••#  4  —  --  4*4  —  *4  ---•••*•  •*•#••--  4*4  *••»---- 

(!•**•,*-•  a. -----  •  4*4  *444— a  Aa  •  »  tt  ••  4  —  —  44  —  4*A*  ••  --a  a  •  »-  ******  ---  *4*  -*••«-  - 

•aflllai-  44  •*••*•••  4  44  4— ■  --*##■•  --  -  444  — ^  44-it-  AAA*  -  -• • fc-  •**•*•  -  -  -  *••  -**•»•- 

.  .ail  4***4*  44**  *44*  4  444  -  -  -  •---•*•*  —  -  44  —  #*.-***•  —  -  4  *»»••  4Aat  A I  444  —  •*-••-- 

..a  *  *•  4  *  *  *  4  44  •  4  — —  4A  A  4  4  -  •—  44  --  44  tt-  AAA  ---  •*  It-  •  A  AAA*  *  *  ----  ••  --  -  a .  4*  - 

•  II  ..  ||||4  ....  •  •  A  A  -•  4  •  A  i  •  -  •—  44  -t-  ♦  >•-•*•••  •—  -4tA-  •**••*••  ---  .  —  -a  •••* 

.4  ....  hi  ..  4  A  A  4  *  -  *44  -**•«•--  —  44  .*#■••  -a  *  a  4* -----  to*  4  A  At  A  A  •  1 1  a  -at*-  - 

t  Al  -a  (a*  44*4  -  4*  A  4  •##**-•••••*•  -  —  -  44  -tt-***---  •  A  A  A  -t--t#t—  4  4  *  AAA  A1*  A  *  *  -  a.’  *  —  — 

. ..  -#t-  4  4  4  AA  *  A*  •  --  A  •  a*  •  -  —  -  ---•*•  —  •**  M  »-•- -a  a  a  -  MlllliM  -it  at*  *-•• 

III  —| . . .  AA  A  4  Mali*  -  ••  — •  *  4A4**  -•  tad  a-  444  •**••*•  -s  tat  -  -• 

.  a.  a. 4*4  --  4*ataA  —  4*4*  - -- *  4  *  tt  44  -  -  4*t*  -*•••*.--  |M»»*M  *(•»-** 

t*a  ieta*  4  Ml  -  a-  4AAt*A  it-  •**•  ---  44-tt  *•  a»M*44*  •  a  a  a  r-  4444AAA*  -••%•  ••  “  a  e 

|i(aa|  ..  4AA4  i  #-  •  4  44  •  tat-IAtaa  -  •  4  aa  -  ••••••  —  a*-  4  —  —  4444AAAA4  -at-  -- 

1  a*  tt  *  4  ....  4  1 A  *  •*  -  a  A  A*  4  a-a  ♦•»*•  ••  4**4  II  •••  -a  44 -  444  *AA\4  ««••—  • 

mu  a  «#t—  444-a-  #•••••  *  *  t  a  •••  --  *AAAa  •••  44  *1 4*  t  —  —  •***•  I  •••-*«-  ta  t  *A 

III  tt*-  4  4  A  A  4  — t  a  4  A  4  44  #  4*AA  •»  A  A*  -  *A  •*•*-•  44*  -4a>  lllllMII  -#*»»--  ••• 

IIM-Itt-  4  *4  4  4t  A  *-•••  -*••»•-  -  “  4  AAA  A-  All-aa  *%•-#-..*,  4A**a4a-  4  aaAAX*  -•••*-  Ml 


•  •••  —  --  Aa  a  4  44  44  4  44  —  4  AA  4  44  --a  a  •  #“ 

•  44  4-a  —  —  —  4  At  A  —  •  444  — ^  44-it-  AAAI  -  -*4t* 

•  a  a  -  -  -  a—  4*t*  —  -  ••  — —  tt'<  -••*•—  *  -*•*— 

•  4-4  —•  4A  A  4  4  -  •—  44  --  44  aa*  AA*  -*-  -t  #t- 

•  «**  -•  4  a  A I  a  -  •—  ••  -t-  *  -•-•*••*  •  —  •*»»• 

•Mil  -  •••  -♦*■••--  —  44  -I* -a  I  --  ••  Aa-»----*  aa  * 

4AA4  4  44  *#--•  at  *A  4  -  — -  44  -tt-***  — -  •***  •••■»♦»  — 

4  *****  •  --  *•«*•  -  mmm  ---•*•  —  *  4  A  A  A  *  -  a- -a  4  *  - 

AAA*  **#•«*  -  ••  —  -•-•  *  4  444*  -a  a*-*  ♦“ 

4*4  --  4  Aa  4a  A  —  ••••  - -- *  4  *  tt  ••  “ -  4*t*  •*••*« .-« 

4  444  -t-  4AA4*A  it-  4AA4  ---  ll«aa  •*  —-•%***•  - a  a  a  r- 


4 a  1  *•••  -  •*•  . . - 

•  4  AA  a  vAAaa  •  a-  .4*  - 


•  III.—  •**•  -•-  ••  I A  A  ♦-•-  4 


tt  4A#AA  aaa*l  - - 

tt-  . . 


Figure  4.6:  Radial  Derivative,  Lower  Quadrant 
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Figure  4.7:  Radial  Derivative  and  Shadows 
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Figure  4.8:  Precursors  of  Shadows 
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MODELING  OF  SHADOWS  IN  RADAR  CLUTTER 


I  Ml  %  SUHt'OUT  I  »4*  WMIIS  r**n*4rrni5,  fr»M‘»UI»n  Gf'lMASC  AN1  C4IIM  MAJOR 
•jK.'lt  'll',  *1  »**’*»  104  A  NIH-Ill*  AL  ClUSI  L  M I  MG  J'JI*  USING 

oni  nr  fur  m<  i.uns  pnhgpamii  •»  is  a  version  »if  guhmuuiihe  "kmi.  an" 

every  jou  iMisuifffs  i igfr  supplied  orcic  uccupnh, 

I  •  program  •wivrn"  performs  i%f  following  tasks* 

A*  ASSIGNS  IN.*Ul/lMI»l».|f  UNITS* 

H.  rGI  AMI  | '.MAS  1  Hf  f»|Nr*|S|ON  OP  f  HP  "  X"  ARRAY  AND  SETS  THIS 

ninths ii'N  rn  "limit* 

e.  talus  succour  inc  *#  x»;C2"  iformcnlv  "CxEC"! 


2*  suwnuusr  "use*"  is  CM»LOvro  to  rcao  the  complete  set  or  scores 

ON  THE  V  All  I  A  PL  fcS  FOR  ONT  OAfA  UN  |  T  * 


1*  ruNCTION  “01  ST-  rnna^Tf*,  IN*  >  |  STANCE  Rt'TWFCN  Two  OATA  UNITS  CM 
OTTW rru  A  OATA  UNIT  AN  >  A  C.UJT*R  CCMTSO|0.  IMF  USr*  CAN  S»EC|FV 
ANY  OFSINCO  01  STANCE  FUNCTION  AHO  WEIGHT  THE  VARIABLES  IN  ANY 
MANNER* 


NOTE  THAT  SCALING  A NO  TP A NT r H9W A T ION  0*  VAO| ARLES  CAN  PR 
ACCOM»LlSMro  EITHER  IN  SUOTOJTINS  "USC9"  09  |N  SUBROUTINE  »OIST«* 


NOTE:  MODIFIED  SURPOlTINE  «*CX*C"  IS  CALLED  "EXEC*"*  ANO  CALLS 

"RSULT2"*  "KMC  AN"  NAS  NOT  SEEN  ALTERED* 


INPUT  SPECIFICATIONS 

CAPO  I  TITLE 

CABO  2  PARAMETER  CARO 

COLS  1-  A  NE-NUNRCf,  «jp  ENTITIES  COATA  UNITS! 

COLS  9  —  1  3  NV*NUMfJER  OF  VARIABLES 
COLS  14-IG  NC*NURRE*  OF  CLUSTERS 
COLS  19-20  NT  IN*  INPUT  UV I T  <0T  THE  OATA  SET 
NT  1 N*S •  CAPO  RFAD59 
NTIN.NF.5.  TAP?  or  DISK  FILE 

COLS  21-2?  NT  OUT*0'J  T  PJT  JN  |  T  <09  SAVING  CLUSTER  MEMBERSHIP  LISTS 
NT  OUT  *7*  CAPO  PUNCH 

NTOUT  *Lr  •  0*  *)D  NOT  SAVE  ME MAER SH| P  LISTS 
COLS  2 3-2?  MINRCL«TFRM INAT I3N  PARAMETER.  CLUSTERING  ENOS  WHEN  A 
CTCLF  THROUGH  THE  OATA  SET  RESULTS  IN  "MINREL" 
OR  EE* ER  CHANGES  IN  CLUSTER 
MINRCL.LC.O.  ITERATE  TO  CO **“*l ETC  CONVERGENCE 
COLS  20—29  |PART*INITI  A.  PARTITION  BA*>\MrT6R 

| PART*! .  SECT  POINTS  ARC  SELECTED  FROM  THE  OATA  UNITS* 
READ  THE  SEQUENCE  NJ4RERS  EOR  THE  CHOSEN  OATA 
UNITS  FOR*  CAUD(S)  .1  IN  20IA  FORMAT.  |F  THE 
OATA  SET  IS  NOT  STORED  IN  CORE*  THE  LIST  OF 
SEOJENCE  NUMBERS  NJST  BE  IN  ASCENDING  ORDER 
IPART *2*  THE  DATA  JNtTS  ARE  GROUPED  INTO  AN  INITIAL 
PARTITION  IN  THE  IN*UT  SEQUENCE  WITH  THE 
FIRST  "NUMB7<1>"  IN  CLUSTER  I*  THE  NEXT 


COLS  30— 31 


"MUNB»|?|"  IN  CLUSTER  2  ETC.  RE AO  THE 
"NUMSR"  ARRAY  FROM  CARD! SI  X  IN  201%  FORMAT* 
|OART*3«  T ME  SCOTE  VECTORS  E DR  THE  SEED  *OINTS  ATE 

READ  FROM  CARO(S)  %  IN  FORMAT  **FHT**  WHICH  IS 
READ  FROM  CARD  3* 

METHOO«*>ARAMf  TEP  E3R  CHOOSING  THE  ALGORITHM  IN  ONE 
VERSION  0T  SUBROUTINE  "KMEAN* 

METHOD*!*  JANCEY  ALGORITHM 
MCTHOO.NC.I.  *  ORGY  ALGORITHM 


•  ♦  CAPOS  3  %N0  A  ARE  Re  AO  IN  SU3PDUT  INC  "KMC AN"  ACCORDING  TD  THj 

•  •  "RODEOURE  SPCCiriEO  HY  ME  CHOSEN  VALUE  OF  "IP4RT".  NOTE  THAT  THE 

•  •  PASIC  K-WEANS  METHOD  OF  w  AC  QUEEN  SI%UY  USES  THE  FIRST  "NC"  OATA 

•  •  UNITS  AS  CLUSTER  SEED  POINTS  ANO  THEREFORE  IGNORES  fMg  **|  » ART* 

*•  PARAMETER* 


STORAGE  ALLOCATIONS  IN  THE  "X"  ATR4V 


X(N|  ) 
X(N2) 

I|K3  I 


TO  X IN?—  I  I 

TO  MN1-I  | 
TO  XfNA-l I 
MM)  TO  XI NS- 1  ) 
XINST  TO  XIN6-I > 
X  (  N6  )  TO  X<  NT | 


NCANV  WOTDS--STOMAGC  OF  THE  CENTR  ARRAY 
NC  WOODS- -STORAGE  OF  THE  NUM.XR  ARRAY 
NF  WJH3S- -STORAGE  OF  T  HE  ME  NOR  ARRAY 
NCtNV  WOBDS — STORAGE  OF  THE  TOTAL  ARRAY 
NV  OR  Nl»MV  WORDS — STORAGE  CF  THE  OATA  ARRAY 
NE  WORDS- -STORAGE  OF  THC  LIST  ARRAY  IN  "RESULT* 


THE  STORAGE  OF  THE  "LI ST"  ARRAY  AFTER  THE  "DATA"  ARRAY  NOW  INO.IES  THAT 
"MAX"  MUST  RE  THE  VALUE  OF  "NT"* 


Figure  A.1:  Forgy-Jancey  Comments 


MODELING  OF  SHADOWS  IN  RADAR  CLUTTER 


CAim  ? 


c 

c  m  g 

1-  R 

c 

C  l>L  S 

9-1  * 

c 

COLS 

1  A—  1  A 

C 

c 

c 

COLS 

I9-.TS 

c 

c 

COLS 

21-32 

C 

c 

c 

c 

1  COLS 

23-27 

c 

COLS 

2A-29 

INIHJI  '.VCIHCAM/WS 
f  MM  t  V  I  Tl  I 

i*Ai*4Mrrn>  r.  apo 

ur  u  or  runilf.  ioataoniig) 

»iv«Nii«*nru  nr  v4i<i  ahles 
O#  CL*/.  IFRS 

MllN-lrP'JI  ij'MI  •  irt  INC  OATA  SET 
Ml  IN-*.,  CAI»M  H‘  A)M 

NiiN.Nr.s*  f**r  tu*  ois*  rite 

Nfrmf*n<n»ur  unit  savi  ns  cluster  mfratrsmip  lists 
NTOt<r>7,  CARD  PtJMCM 

NtW#T  ,|  r on  »mf  SAVE  membership  LISTS 

M|NM* L*T<?UW|N4r  I  <)H  PIBAir-’ffW,  CLUSTERING  ENOS  WHEN  4 
CVCl.  F  fMNOiKiN  THE  0*14  S^f  wrSULIS  IN  -M| NREL* 

OH  Ft#FR  CHANGES  IN  CLUSTER  Mr.  HOC  ASM  I  PA 
MINMCL.Lr#0,  ITrpATr  TO  fOW'ir  TF  CONVERGENCE 
IO  APT  *  !Nf  T  I  A.  >ADffTf<V«  ^AQA^rffA 

|  PARI  *  I  ♦  POINTS  Apr  sr^rcfro  FROM  THE  DATA  UN  ITS. 

PE  AO  TMr  SEOUCNCC  NJM9CRS  EON  THE  CHOSEN  OATA 
UNITS  FORM  CAPOISI  **  IN  20  14  POP^AT,  IF  THE 
DATA  SET  (A  NOT  StflOFO  In  COPE,  TMC  LIST  OE 
SrojEMTE  NUNOEPS  -JST  fl*7  IN  ASCCNDlNG  ORDER 
I  PART  *2*  TMr  OATA  UNITS  APE  GRCU*FO  INTO  AN  INITIAL. 

PARTITION  IN  in'*  INPUT  SEQUENCE  WITH  TMC 
FIRST  -HUM«**«|»-  |rl  CLUSTER  |,  THE  NEAT 
-NUW0»I3|-  IN  CLUSTER  2  FTC*  REA0  THE 
"taUMQtt*  AMPAY  FROM  CAPOISI  3  l»l  >014  FORMAT* 

TMF  SCONC  VECTORS  FOR  THE  SCED  POINTS  ARC 
REAS  FROM  CAPOISI  4|N  FORMAT  -FHT-  WHICH  IS 
RFA3  TROW  CARO 

MCTHOO*PARAmb tpn  FOR  CMOOSINS  1>E  ALGORITHM  IN  ONE 
VERSION  O*  SUBROUTINE  -XMC  AN¬ 
NEXED*!  •  JANCEV  ALGORITHM 
METHOD. nE.I *  «ORSV  ALGORITHM 

NCWAX««UX  NUN RF W  OE  CLUSTERS  FOR  OUTLIER  ROUTINE 
NCNAXvNC*  DON«T  00  OjTLlEP  ROUTINE 

NC MAX  *LE .  0.  UST  NE  AS  MCM4X  fl.F*  NO  LIMIT  ON  CLUSTERS! 
OUTL VRrCOEEF I C I ENT  IOFTWFFN  )  AND  t  INCLJSlVEI  FOR  THE 
OUTLIER  ROUTINE  (NUMRCR  IN  F4.2  FORMAT  I 
OUTLVRbI**  MAX  AMOUNT  OE  OUTLYING  CLUSTERS 
OUTLVR****  NO  outlying  CLUSTERS  ALLOWED* 


I  PART  *3, 


C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

CPPV  CAROS  3  ANO  A  ARE  RCAC  IN  SUOPOUT INC  -*ME4N-  ACCORDING  TO  THE 
CT**  PROOEOURE  SPEC  I  El  CO  MY  TM«  CHOSEN  VALUE  0-  -|P4RT".  NOTE  THAT  THE 
€•••  I* ASIC  W- MEANS  MFTHCO  OF  MACOUF.SH  S1M*»LY  USES  TMT.  FIRST  "NC*  OATA 
C*#p  UNITS  AS  CLUSTER  SEED  POINTS  ANO  THEREFORE  IGNORES  THE  "IPART* 

C#P#  PAMAMETER* 

c 
C 
C 
C 
C 
C 

c 
c 
c 
c 


COLS  30-31 


COLS  32-39 


COLS  AO— A3 


STORAGE  ALLOCATIONS  IN  THE  "X-  ARRAY 


X|M  I  TO  XIN2-I  I 
XCN2I  TO  XIN3-I I 
XIH3I  TO  XIN4-I ) 
XI HA |  TO  XI NS- 1  I 
XtNSI  TO  XIN6-I I 
XCN6!  TO  XIN7I 


NCPNV  WO>OS — STORAGE  OE  THE  CENTR  ARRAY 
NC  WORDS-— S TO R4G1  OF  THE  NUM3R  ARRAY 
NE  WORDS- -STORAGE  0-  THE  ME*RR  ARPAY 
NCPNV  WC»DS — STORAGE  OF  THE  TOTAL  ARRAY 
NV  OR  NE4NV  WTRDS--STORASE  O*  THE  OATA  ARRAY 
NE  WORDS- -STORAGE  OF  THE  LIST  ARRAY  IN  -RESULT* 


THE  STORAGE  OF  THE  -LIST-  ARRAY  AFTER  THE  -OATA-  ARRAY  NOW  IMPLIES  THAT 


£  •MAX"  MUST  OE  THE  VALUE  OF  -N7-* 

CPPPPtPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPP* 

c#  . 

c#  note:  1 

c#  since  nc.  can  ne  incpcasso  py  the  outlier  routine*  the 

CP  ALLOCATIONS  IN  THE  AMRAYS  ARF  CASED  ON  NC  MAX  *  NOT  NC* 

c#  THIS  CAN  INCREASE  ST0RAS2  REQUIREMENTS  FASTER  THAN  EXPECTED* 

C* 


'Figure  A. 2:  Outlier  Comments 
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FORTRAN  IV  G  LEVEL  21 


exec 


OATC  •  T934S  09/54/50 


0001 


0002 

0003 

0004 

0005 

0004 

oo  or 
0000 
0004 
0010 
001 1 

0012 
0013 
0014 
0015 
0014 
001 1 

ooto 

0010 

0020 

0021 

0022 

0023 

002  4 

0025 


•SUBROUTINE  EXEC  I X •LIMIT} 

C  THIS  SUBROUTINE  REARS  PARAMETERS.  COMPUTES  STORAGE  A NR  CALL'S  MAJOR 
C  PROGRAM  SEGMENTS  NF*DE9  F0»  A  NQN-M I FQ  ARCH  |  AL  CLUSTERING  JOS  USING 
C  CNE  OF  THE  MgT MOOS  ’ROGRAMMEO  AS  A  VERSION  OF  SUBROUTINE  ■KMEAN" 

c 

C  MORE  COMMENTS  AND  A  FULL  EXPLANATION  CAN  QE  FOUNO  IN  FILE  •OO.FCCMMENT 
C  OR  nv  OBTAINING  THE  REPORT  THAT  MAS  A  SOURCE  OF  THIS  PROGRAM. 

C 

DIMENSION  X|  |  I.TITLEI20I 
READ  15.10001  TITLE 

RC AO  (S.ltOO)  NE. NV.NC .NT IN.MTOUT  » MlNREL . I PART .METHOD 
•(RITE  16.2000)  TITLE 

WRITE  16.2193)  NE  .NV.NC.NTIN.NTOUT.MINRCL. IPAR  T. METHOO 
Nl*l 

N2«N| 4NCANV 

N3«N24NC 

N4sN34NE 

NS«N44NC4NV 

C  «N6*  MAY  CANO  WILL)  BE  INCREASED  IN  •UNBAN" 

N6»N54NV- | 

N7-N44NE-I 

MAX«N6 

WRITE  (6.2200)  MAX, LIMIT 
I F  (MAX  .GT.  LIMIT)  STO® 

CALL  KMEANC  Xf  N|  ).  X(N2|  .X  (  N3)  •  Xf  N4  )  .  XI  NS)  .  NS.  N€  .NV  ,NC  .NT  I N.  M|  NR  EL. 

•  |»ART.M£TMOO.L  IN  IT ) 

CALL  RESULT  I  X (Nl ) • X (N2 )• X ( K3) ,X(N4 ). T I TLE • NE.NV. NC.NTOUT. 

•  X(N4).X(N5I) 

METURN 

1000  F0RMATC2944) 

1100  FORM  ATI  SI  S ) 

2030  ECRMATH-M  .29  A%  ) 

2100  FORMATCONF  *M1,/.»  NV  **,|(t./,«  NC  a*.  !»./.•  NT  I N  ••,16./. 

•  •  NTOJT  -».IS./.*  MlNREL  -*.|4./.*  I® ART  •‘.IS./.*  MCTMOO  »*.I4) 
2200  F CRM AT ( * 3RE0U IREO  STORAGE  »•.)»,»  WORDS*,/. 

•  *0 ALLOTTED  STORAGE  •••IS,*  WORDS*) 

END 


Figure  A. 3:  Subroutine  EXEC 
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FOOT  RAN  IV  G  LEVFL  *1 


MAIN 


OATS  »  79349 


0002 

0003 

0004 

00  OS 

0006 

0007 

0000 

0009 

0010 

001  t 

C  "1 

0012 

0013 

0014 

0015 

0016 

0017 

0010 

0019 

0020 

1000 

0021 

1100 

0022 

2000 

0023 

2100 

0024 

2200 

0029 

ALTERED  SjnROJTINS  EXEC — CALLS  KMEAM#  RSULT2 
SUHROUTINE  EXFC2IX, LIMIT) 

THIS  SUUftOUT  INE  RE  ADS  P4B*  4E  f  CPS  •  COMPJTFS  STORAGE  AND  CALL •  S  MAJOR 
PRCGRAM  SEGMENTS  NESD*D  FOR  A  NON-Ml  F«A»CH|  AL  CLUSTERIN';  JO?  USING 
ONE  OF  THE  METHODS  PROGRAMMED  AS  A  VERSION  DF  SUQR3UTI ME  "KNEAN** 

MOPE  COMMENTS  AND  A  Flj__  EX  >L  ANAT I RN  CAM  «?  FOUND  IN  FILE  ••DO.FCCMWENT"# 
OR  PV  OBTAINING  THE  REPORT  THAT  WAS  A  SDJRCS  OF  THIS  PA OS4  A 4 • 

— CHANGES  MADF  ID  -N#,»,  SURPOUTIMF  -RESULT*  IS  "MSULT2*# 

- — CHANGED  IM^UT  * OPM ATS  (DOCUMENTED  IM  ROSCDE  FILE  "FCOMMENT**) 

DIMENSION  X(  »  >.  Tf  TCEf  20) 

READ  (StlOOO)  TITLE 

HEAD  I S  #  I  I  03  I  NE.NV.NC #NT| N.NTOUT •MINREL .  IPART.NETHOD 
WRIT*  (6.200?)  TITLE 

WRITE  (6,2100)  NE«NV.NC»NTIN.NTOOT«M|NREL»|PART*METHOD 
N»»l 

N2-NI ♦NC*NV 
N3«N2+NC 
N4«N34NF 
N5«N4 ♦NCANV 

46"  IS  NOT  THE  SAMS  AS  IN  ORIGINAL  PROGRAM#  ALSO  ”N7"  CHANGED#  TOO 
N6*N5 ♦NVANE 
N  7  *  N6  *  NC  -  I 
M A X*N7 

WRITE  (6# 2200)  MAX, LIMIT 
IF  (MAX  . GT •  LI  MI T |  STO» 

CALL  KMEANI  X(  Nl  )  •  X(  N2  >  •  XI  N3 1  •  XC  N4  )  #  X(  NS  )  •  NS«  NE  *NV  #  NC  •  NT  IN#  M|  NREL# 

6  |PART,  M£THOD#.  IN  IT) 

CALL  RSU.T2(«(Nt) #X(N2)#X(N3) #X (N6) #T ITLE.NE.NV# NC.NTOJT# 

•  XINA )#X(  NSI) 

RETURN 
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NRrPaOtrtf 
JARPtJ 


Figure  A. 7  (Part  2):  Subroutine  KMEAN2 
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I'MtHAN  IV  ft  LCVri.  2 | 


ftftTAN* 


Mvr  •  mu 


•vm/ii 


soao 

00*1 


6062 


00*13 

60*4 

00*6 

00ft? 


140  fCNriwf 

A VII  *il  ■  fWNtM/IIC 

c  test  rn»  ivuur's*  and  auv.atc  ntv  cufftTron  I*  allow* 
ir  t  .n*h«  *LP*  Avoiniiniir.vHi 

•  .  ,tfO,  (NC  *LV«  OCtMIII  60  TO  |4ft 

C  NFOJ  NTtt  CCtffttrM 
Nr*>«c»i 
JOFFINC 
r»r  14ft  L«|*NV 

(rNnMNVHiRrr-ii»il<OArMK|»LI 

CflNf  IMiF 


141 


0066 

1004 

0000 


0001 

0041 

004) 

0004 

001ft 

0006 

0097 

0096 


IftO 


0090 
0100 
9101 
0102 
0103 
.  0104 


010ft 

0196 

0107 

OIOO 

0100 

0110 

out 

0112 

0113 


0114 


0119 

0116 

0117 

0110 


0110 

0120 

0121 

0122 


AlLOf  AVC  OAT 6  0*1  If  •«"  TO  CLUSTF6  "OOFF* 
l<U«"tW|  J«"F  I.WJNfint  JH*F  141 

I A  fj  urp  #r  o.  *r*naf«||  40  TO  IftO 
IMF  OAT  4  UNIT  CMANftgft  ITS  NEN6FF  SM|» 

A0VT4»4WE5M 
NF  4***f«l*  J*Cf 
Jl*l JR€F-||4NV 
OO  160  l«I.NV 

«t*«|4| 

TOTAL  I Jll •TOTAL ( J| |*9ATA|K| » 

160  CONTINUE 

C  ALL  OAT  A  UNITS  ALLOCATE*  TEST  006  CONVEOOCMCC 
»PfT?|6«*700>  M0vFS*N7A)4«f0)Sf*NC 
N»AftS«N»ASSAl 
J»CF*0 

IF  (MOVES  *6T,  6IN6CLI  ftfl  TO  Iftft 

IF  (METH3D  .NC.  |  .ANO.  MOVES  *10*  01  MflNN 

JPEF-I 

C  COMPUTE  TRU*  CL 'A  ST  £4  CCNT6OI0S~FO6CV  U»MTI 
170  Jl«0 

00  ISO  JM.NC 
90  100  1»I.NV 

IF  INUMD4IJI  .Lf .0)  CO  TO  160 
CFNTrtf J| |»TOT A.(j| l/NURORI Jl 
too  CrKTlNUC 

ir  IJRCF  .CO.  II  RETURN 
60  TO  129 

109  IF  INETH03  .NT#  1 1  GO  TO  170 
C  JANCFV  UFO  AT  C 
190  J1*0 

DO  700  Jtl.NC 
00  200  I* I *NV 

C  TEST  FOP  NUM&U0  ►«**•  TOO... 

203  CENT6IJI l*2.AT0TAL( Jl I/NJNS6I J1-CEN?6| J| » 

60  TO  120 
FC6VATI29A4I 
FC6NATI20I4I 


1090 

1100 


F06T6AN  IV  ft  LEVEL  21 


R6CAN2 


DATE  «  2036ft 


00/ft 4/ftl 


0123 

0124 


0120 

0126 

0127 

0129 

0120 

0130 


0191 


2000  trPW3^m5°6?4CLUSTE4  ANALYSIS.  OATA  KT  ST06EO  IN  C06E»I 
2100  F06NATI  *3PF0UI6E9  STORAS*  «MS**  NORMS/* 

6  *0 ALLOTTED  ST06A6E  «SI1*«  60609*1 

FOftMATt  *0/1466  AT  •  «  26  A4  I 

FDTVATf  * I  INITIAL  C LUST £6  CENT  EPS  RCAO  IN  AS  EOLLOOft*  *///! 
f 004 ATI IN .1 6« 12.41 

FC6MATI «|  IPAPT  «••(!•%  NU«OR  A66AV  6240  AS  FOLLO»ft*.///l 
enPMATC  IX .  I  01 71 

FOWMATI |N0. IS.*  DATA  UNITS  M9VC0  ON  ITE6ATI0N  N/NAEt * . I  3. /• 

•  •  SUMNO  DEVIATIONS  AftOUT  SEED  POINTS  « • • Et 6. »• I OX. «NC  »*.ISI 
ENO 


2290 

2390 

2490 

2900 

2600 

2700 


Figure  A. 7  (Part  3):  Subroutine  KMEAN2 
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fwmw  iv  «  irvn  1 1 


MIR  •  HIM 


MiftiHw? ivr  rr sult icf nv b* mi vtf** nFvim#i, ist.t itlc **r.Nv*Ne*NTOur* 

A  fill  A4  *'»4«A) 

c  ihii  wtmnjt  ixr  inf  pcwh  tbon  a  ausiHm  jm  basis 

€  ON  ANY  VfcWMON  rtf  ViHmuiYI Nit  "KMC  AN". 

€  CCNf Ml  1 1  *NJNB»I  I  I  aMEMNR|  ll*L ISM  lMIVlC(tB|»VOTAMIIi 

•  DATA! I I 

C  AS  A  CfMflMVNCI  CAIT  (ON  ROUE  OUT  THC  BAV  NCMilSMSMfP  LIST. 

TITL* 

iPlirlA.MB^I  |ntn<%N«r»»ic«I9NI!I 

vniiriAfim)  f  MjwH'ti  J  |  •  Ja  |  *NC  I 

C  INV»  RT  I  MR  "BCMK-  ARRAY  A NO  »Uf  f  MV  BfSUlT  IN  THE  %IST«  A«BAV. 

c  riRsr  orvisr  ?»r  a^rat  io  contain  iiaii  points  in  the 

C  "U'.T"  AHflAV  VON  CAT.M  CLUSTER 
NUMON INCI  «NE-  MUNIS  INC  I  VI 
JJ"NC 

no  19  J«>|NC 

NUMMlJil  MMIMUMJ^-NUnKiill 
JJ-JVI 
JJIWJJ-t 
|0  CUNIINU* 

C  BUILD  "LIST-  ARRAY 
OO  79  *  — I  »NV 
ME  MONA*  NCNORIKI 
NjatwMOtlBCNQMl 
LlSfCNJIaK  ^ 

NUNOAI NVNBNKI  vNUNBUlNCNBim  I VI 
10  CONTIMJV 

C  SAVr  THE  snurro  NCMOCNSMIP  LIST  IP  OCSIRCO 
IV  INTOUT  •l'«  91  M  TO  30 
WRITE  INTOUT*  39991  TITLE 
WRITE  C  NT  OUT  *  ?|  09  I  «L  I  STI  K|  *Kw|  »NC» 

C  RESTOPC  THE  "NO-ABB"  ARRAY* 

39  UJ"NC 

OO  A9  J"«*NC 

NUMB  I JJ »  MUBBAUJMUNBBI  ii-ll 
I 

AO  CONTINUE 

NUNN* (Its NUN**! ||*| 

C  PRINT  RESULTS  FOR  EACH  CLUSTER 
WRITE! A «t 000)  TITLE 
Kiwi 

OO  SO  J"l  *NC 

WRITE! 6* 9 300)  J.NUNBNI Jt 

Jl*« J-l  )•  Nv 

WRITE (6.E A00)  ( CENTRf  J|V| 9  *l»l *NV) 

KtwKI VNU3 SR I Jl- » 

WRfTeC««?S009  TLf Sf  TK I »K*K I *K£f 
KIwKSVl 

?°  TMll°SHpSeOUT  "TOTAL"  ARRAY 


FORTRAN  IV  «  LEVEL  SI 


DATE  •  HIM 


S?,f?«JTli8S>  i 

VORWATI*-TOTALS  "OR  CLUSTER** I 3) 
«RITCI6*AIOOMTOTAL<tVJ-|)  *J»I»NV) 

FORMAT! IX.EtO.S) 

CONTINUE 

RETURN 

VORWATf  |H|*E9AA  I 

V CRN ATI *9 RAW  ME WAEBSNIB  L I  ST* ./ . 1 1 X.fS IS M 
rORNATI  *0 CLUSTER  IU*S*«/»nx*HISM 
FORMAT! •^CLUSTER* • 1 3* *  CONTAINS** IS* *  DATA  UNITS* I 
ropMUHCTNTwn  Coordinates* */*(tx«tocis*Ai) 

VORNATI •ONCPAEASHI"  Ll ST*  */•( IX.2SIS)  » 

FORMAT! SO A4) 

VORPATIIBIS) 

CNO 


Figure  A»8:  Subroutine  RESULT 
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UHlMAtt  IV  •  mVFk  >1 


Mtt.ll 


mu  •  mu 


t*/}|SIS 


0001 

c 

c 

c 

0009 

0001 

0004 

0000 

c 

c 

c 

c 

<1000 

0007 

0000 

0000 

0010 

0011 

0019 

0019 

10 

c 

0014 

001 S 

0010 

0017 

0010 

OOIO 

90 

c 

0090 

0091 

0099 

002  J 
0094 

99 

c 

0090 

0090 

90 

0097 

0090 

0090 

40 

0090 

c 

0091 

0099 

0099 

0090 

nn 

••37 

••;• 

••j* 

I  ill  M.  •  IH  I  A  I 


THIS  MAlWMiriV  f Ml?  NFAm  in  IRON  A  CLUtURlNO  /W  IA1ID 

on  mu  vcmsion  if  sjmrijvflnr  h«i:rn«« 


DINFNmLn  »NTR(||(N*INM||)iMrM|M||)iLlinMtVIVlMtlMOrAUlli 
•  MfAIII 


An  a.  CONflN'.fHC*  AMC/MiriON  «ft|Tt  mil  IMS  Ml  NfMASItNl#  Lift* 

»p|T»*in*?ooo>  title 
«<>ITf  <0*21001  l«r<n»(Rh«i|»liei 
(NI^OhJvliNCI 

INVi  Pf  I K  4IIUAV  AND  MM  INC  BfSULT  IN  INF  •LIST*  A«BA 

l|»F,f  AC VI nr  TM?  •VJVMP  AONAV  TO  CONTAIN  STABT  POINTS  IN  THE 
ANN  AT  FO*  f  ACM  CLUSTF.B 
NO NORI NCI  *NE-KJN'JB|NCI  *| 

JJ«NC 

DO  10 

NUVIMU4I  IsNOMIOl  JJt-NUNBBf  JJI) 

JJ9 Jit 

CONTINUE 

BUILD  -LIST**  AAA  AT 
no  2»  Ml  *NC 

NCN3BK«M?NftB| K| 

NJtW|MM(«WIAO 

niriNiiti 

NUNAA I NCVBfNO  •NMMIRENMK  |»| 

CONTI  Nur 

SAVE  THE  somo  MEN1CASHIP  LIST  IF  DC  SI  BCD 
IF  INTOUr  .If*  1|  <50  TO  SO 
fcAlTE  I  NT OUT *39 00 )  Tt  TLP 
•  BITE  | NT  OUT • 31 0) |  (Ll ST< K» .«■! .NCI 
DO  22  R«liNE 

CALL  USBUUTOATAIILISTIIO-IIANVAIU 
BCSTOBE  THE  •NUNBB"  ABBA V  * 

Ji*NC 

DO  AO  i»2«NC 

NUniMi2)iNU<iaR(jjMU«RU2-|| 

JUOJ-I 

continue 

Nil  NAB  (  I  I »  NUNBBI  1 1  - 1 
POINT  NESULTO  FOB  EACH  CLUSTCB 
«B|TEI0*7000t  TITLE 
Rl«l 

00  SO  i*l  .NC 
VBITEIA.FSOOI  JiWINMUl 
JI«(i*IIMB 

»PfT£ 16*7400)  f  CENTO! 91*1 1  *!•! *NV| 

K9«Bt»NU«8B(Jt-| 

VBITEI6.2500I  <LtST|K)*K«Kl*K2t 
KKR2N 


FOOT  BAN  IV  0  LEVEL  tl 


BOULT t 


DATE  ■  70901 


t 0/91/1 9 


0040 

2*  T*M°SHS,S(»T  wmt 

0041 

90  00  1-ltNC 

0049 

VBtfPIO *40001  1 

0049 

4000 

FOBNATt •-TOTALS  FOB  CLUSTEB**I9» 

0040 

»B|TE(0*4 1001 1  TOTAL!  U-t  I4NVF9I  •  JM  eNVI 

004  S 

4100 

F0BN4TtC|K.Et0*S»t 

0040 

00 

CONTINUE 

•047 

BE  TUB  N 

0940 

9090 

FOBNATI  Ml  *70 A4  | 

0040 

9100 

FOBNATI *0B4V  NSmOEBONIP  Ll ST* •/•( 1 K *25151  9 

0090 

9700 

FOBNATI *9CLUSTEB  SITES**/*! IV .751511 

0051 

2900 

FOBNATI •OCLUSTCO*  *  1  ?*•  CONTAINS* • IS* •  OAT  A  UNITS*  I 

0059 

2400 

FOBNATI*OCENTB9|.n  C09N0I NATES*  */*HR*  1 0019*41 1 

0059 

9509 

rOBNAT(*9NENlCBSN|P  LIST**/*! IX *951911 

0054 

9099 

FOBNATI  72  A4I 

•055 

9100 

FOBNATI  15 151 

•050 

ENO 

Figure  A. 9:  Subroutine  RSULT2 
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I*  a  LLVCL  II  4/  IS  MIC  •  ••III  -  ll/lt/M 

C  o«iv*r.  rutii  tci<  use*  in  filI  •fcunmsmt* 


0091 

OUCNSIO.  Nl  14300 1 

0091 

LIMIT-24900 

9903 

CALL  tXEvJM.LlNU) 

0099 

*N*>  6  |LC  4 

0003 

S»TO*» 

0006 

CNO 

FORTRAN  IV  9  LEVEL  II  OUT 


0001 

FUNCTION  atSTfX'V) 

0002 

9mN5in:  xuiivni 

c 

DISTANCE  IS  EUCLIDIAN  fX-Vl 

09)3 

OlftT-O.9 

0004 

DO  1)  l-|  ,< 

000ft 

to 

91  ST -01  STAC XI 1) -71 1» |462 

00)6 

AC TUMN 

0007 

ft  NO 

•41*  «  Milt  04/59/04 


r3HI04N  IV  «  um  *1  “•!•»  "" .  RATE  •  40114  04/41/10 

3TTI  EVNCTIO-)  OISTIX.VI 

eui  9i»i.45iC'i  xiii. rm 

OOTJ  C  ’"^lIKftiiirSlSlSfift  C,*tV  “***  THe  ""**  VAmAW-INTINUtT 

00  5%  *  iH’TU4A 

CJOft . .CNO 


Figure  An  11:  MAIN:  Two  Versions  of  DIST 
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FORTRAN  IV  C  LEVEL  « I 


OATS  •  7934ft 


DIMENSION  VALUE  (  |  00*  100) 

S I GMA |S|« 

S I6M42* 10000. 

I  X*1 23459  TOO 
00  100  1*1.100 
no  IOQ  3*1.100 
CALL  RAVI  IX.lV.VFL.StONAtl 
100  VALUE I I »  J) *YFL 

on  I  1*1.  IN 
no  i  J-1.100 

CALL  R4VMX.|V.VFLtSI6Nllt 

1  VALUE! f.Jl* VALUE!  1.31* VFL 
OO  ?  1-73,30 

00  2  3*1.39 

CALL  RAVI !X.|V.VFL.SI0MA2) 

2  VALUE  I  1  •  Jl -VALUE! I  .3)  + VFL 
00  3  1*23 • 30 

00  3  3*M  .100 

CALL  RAVI  IX.l  V.VEL.SI6MA2I 

3  VALUE!  I.JI-VAL3E!  I.3I+YFL 
OO  4  |*3I *99 

no  A  3*1. 100 

CALL  RAV|  IX. IV.VFL.SI0MA2I 

4  VALUE  1 1  •  -I)  -VALUE!  I  «  31 A  VFL 
00  9  1*60.40 

DO  9  3*1.99 

CALL  RAVI  IX.I  Y.YFL.SIGM42  » 

9  VALUE! f.3)*VALUei  I.JIAVFL 

0*  a  |*60 .60 
DO  6  3*01 .100 
CALL  RAVI tX.IV.VFL.5IGMA2» 

6  VALUE! I .3) -VALUE! 1 .31 A VFL 
OO  7  |*6| .100 

00  7  3*1.100 

CALL  RAVI  IX. 1  V.  VFL.SIGMA2! 

7  VALUE! 1.3 I -V4LJE! 1.31 AVFL 
OO  300  1*1.100 

00  300  3*1.100 

309  WRITE! 1.9001  I .3 .VALUE! V* 31 

900  FORMAT 1214. E 17.1 01 

END  FILE  1 
VRITEI6.I0001  IV 
1000  FORMAT! IX . 1 1 1 1 
ENO 


Figure  A, 13:  Creates  Rayleigh  Data,  SIGMAslO.OOO* 
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FORTRAM  IV  «  LEVEL 


OOOft 

000* 

0033 

0004 

0003 

0006 

0007 

0000 

0009 

0010 

00ft  ft 

100 

001* 

0013 

0014 

0015 

1 

0016 

0017 

0010 

0019 

2 

00*0 

00*1 

302* 

0023 

3 

00*4 

00*5 

0026 

0027 

4 

0020 

0029 

0020 

0031 

5 

003* 

0023 

0024 

0039 

6 

0036 

0037 

0020 

7 

0039 

300 

0040 

004  1 

004* 

300 

004  3 

0044 

1000 

at 


MAIM 


01  I  ON  VALUE  f I  00*  I  001 

ftftC'tAIMt 

SIC*442«S9. 

t  **40224/049 

OO  100  1*1.100 

00  100  Jsl'lOO 

CALL  RAVI  | X.  I  V.  VFL .51  <#NA |  ft 

VALUE ft  I.Jft*VFL 

OO  I  1*1.10 

OO  I  J*l.  100 

CALL  RAV| |X « 1 V.VFL.SI 0MA2I 
VALUE! I . J> -VALUE! I »  J) A  VFL 
00  a  1*20.30 

oo  a  j*i.30 

CALL  RAVI IX.IV. VFL.SICMA2) 
VALUE!  I.  J)*VALUEI  l.JftAVFL 
00  3  1*20.30 
OO  3  J*5i  .tOO 

CALL  RAVI  |X.|V.VFL.Sf6MA*ft 
VALUE!  1  .  UPVALUE!  I  •  JftAVFL 
00  A  |»3|  .09 
OO  4  U«l.tOO 

CALL  RAVI  IX.IV.VEL.StOMAal 
VALJF I  I.  JMVALJfC  I  .J14VFL 
OQ  3  1*63.00 
DO  5  J* I . 59 

CALL  RAVI  IX.IV.VFL.SIC«4*ft 
VAIUEIf  .  i)«  VALUE  I  I.JI4VFL 
uU  t»  1*69.90 
00  6  3*91  .100 
CALL  RAVI  IX.  IV.  WL.StCMAaft 
VALUE  I I.J»*VALJCI f.JI »VFL 
00  7  1*81 .100 
00  7  J*ft» 100 

CALL  RAVI |X.tV.VFL.SIO*A*ft 
VALUE!  I  .JMVALUPI  I.UI4VFL 
00  300  1*1.100 
OO  300  U« I .100 
fcRITf(I.SOO)  I.J.VALUCII. J| 
VORMATUI4.Eft7.IOl 
CWO  rftif  ft 
VRITEI6.I  OOOft  IT 
FORMAT! IX. Ill  ft 
ENO 


oate  *  r«m 


Figure  A. 14:  Creates  Rayleigh  Data,  SIGMAs50 
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POMIIAM  IV  0  LEVEL  21 


•LOT 


Oort  ■  rm» 


1324 

0003 

0004 

OOOS 

0000 

0007 

0000 

0000 

1 

0010 

001 1 

0012 

20 

0013 

0014 

0010 

0014 

40 

001  f 

30 

0010 

10 

0019 

14 

0020 

09 

0021 

IS 

0022 

100 

0023 

0024 

0029 

200 

fttHE  S?SSamm<S&lMisoi  .oof 

OAT  A  C*1A»5/*  I  ••  |»J»  •  •0»,«6*,*r#,»0«*»*».»0*.iA».*3*< 

DATA  WT/».«/ 

DO  |  1*1,  IRO* 

DO  I  J*  1  •  JCOC 

aoxc  i,jmoot 

DO  10  1*1  *9999 

READ  (S,?9.*iOt|S)KOU*lT 

FORMAT! IS » 

00  31  NCHOERftlt  (COUNT 

R(!49(  I«40,*N9*I0I  OATA1  •  94TA2 

F0R*4T(»rt.9> 

I NTI • INTI  9ATAII 
I NT?* INTI  04T42) 

00X1  1  NT  1  , 1NT2 1*CNARSC f I 

CONTINUE 

CONTINUE 

WRITER, SSI 

FORNATI*  Ml  ERROR  ill  END  OF  OATA  SET  ON  UNIT  01. •  I 
00  199  9.1.IR04 

«RlTE(10,200M90X<t,Jt  ,J«t  •  JCQLI 

F0PM4TI IX ,12941  I 

STOP 

CNO 


Figure  A. 15: 


Cluster  Plot  Program 
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FORTRAN  IV  6  LEVEL  it 


*A  IN 


OATS  ■  (Mil 


C 

000 1 

0302 

0003 

C  Mi 

300A 

0005 

5 

C 

C 

C  I 

1036 

0-307 

0005 

>030 

0010 

.  to 

001  1 

C  I 

0012 

0313 

001  A 

0015 

3116 

001  7 

0015 

0010 

0020 

0)21 

AO 

9022 

0073 

30 

UFA 

0075 

20 

0026 

0027 

100 

0076 

101 

0029 

102 

00  JO 

201 

0031 

200 

0  332 

A)0 

0033 

300 

00  SA 

644V  SCAL*  *R36*AM#«# 

INtFStH*-'  O^AV 
r  tAL  XL  I'lEC  12.11  •TIWESHO) 

icaiot*i  sc4ie(  n.»u»Eii2i) 

C4i>  IN  T ML  1MAV  SCALE  nVNOOcS*** 

»CAl(5»5l  I  SCALE  (11*1*1  iAI 

fcjrmatioaii 


RE AO  THE  Lf.«EST  VALUE  FOR  EACH  CRAY  LEVEL* 

JCl  I)  l>2tl 

REACH  5.  10  01  THRESHl  1 1 

I  !■!-» 

*P!T«:<A,l)n  SCALE!  It  I  »f  HA  E  SHI  1 1 
CONTINUE 

MM  HE  (6*1  02)  SCAc'il  31  •THRESH!  01 
HCA3  MAX  HJ4|*«  4*0  Vfc*T.  AXIS  VALUES*  «>  lt«t  <-  1  SSI 
WEAO  <5«?30l  *4KV*N44M 
maxv.maxh 
00  20  I  *1  •  MAX V 
DO  30  J»t •MAKH 
RE  AO  11  *3001  X.IMCIJ) 

6*  AVI 
DC  AO  K«M 

IF  | XLl N? ( J)  »LC«  THRESHIRII  60  TO  4> 

6R4VK 
CCMIMC 

t*l  |N?(4)«SCALE(6RAVI 

CUNT1NJE 

*RITC|o*%  00 1 1 RL INC! JJ •  0»l •NASH) 

CCNTINUF 
S  TO* 

FORMATIFIS.T) 

Ft-RM ATI  •  SVM30.**  ••Al.  •••  <««F|S*?» 
r  URMATI *  SVM10L*«*.A| •’••«>"• F I  9 • T I 
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FIGURE  A. 16:  Grey  Scale  Program 
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MODELING  OF  SHADOWS  IN  RADAR  CLUTTER 
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